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Q\ ABSTRACT 

These lectures provide an introductory review of big bang cosmology. I 
discuss the expanding Friedmann-Robertson- Walker universe, summarizing 
the observational evidence which has led to its adoption as the 'standard' 
cosmological model and reviewing its basic properties. Subsequent lectures 
provide an overview of the early universe. The final lectures give an intro- 
duction to the inflationary universe, beginning with the motivating puzzles 
of the standard cosmology (the horizon and flatness problems) and end- 
Q_i ing with the inflationary production of quantum field fluctuations and their 

q possible role in seeding the large-scale structure of the Universe. 

1. Introduction 



In the late 1970's and early 1980's, theoretical cosmology underwent a renais- 
sance: extrapolating concepts from particle physics, in particular the standard elec- 
troweak gauge theory, to very high energies, a framework emerged in which one 
could meaningfully speculate about the evolution of the very early universe. This 
marriage of particle physics and cosmology led to a number of remarkable develop- 
ments, including models for the generation of the baryon asymmetry, the inflation- 
ary scenario, the notion that topological defects could be created in cosmological 
phase transitions, and predictions for non-baryonic particle dark matter, to name 
just a few. 

In recent years, observational cosmology has been undergoing its own rebirth. 
There has been an explosion of information on the large-scale clustering of galaxies 

*Based upon lectures given at TASI '93 (Boulder, CO, June 1993) and at Baksan International 
School of Particles and Cosmology (Baksan, Kabardino-Bulkaria, Russia, April 1993) 



1 



from redshift, peculiar velocity, and photometric surveys gathered by ground-based 
telescopes (see the lectures by Kirshner in this volume). Studies of rich clusters of 
galaxies via their gravitational lensing effects as well as X-ray emission from hot 
intracluster gas have started to provide new clues to the distribution of dark matter. 
In addition, the recent detection of large-angle anisotropics in the cosmic microwave 
background radiation by the COBE satellite provides the first probe of structure 
on very large scales, while a series of anisotropy experiments on smaller scales now 
have tantalizing results. On the scale of the universe itself, there has been steady 
progress in attempts to measure the cosmological parameters (in particular, the age, 
expansion rate, and mean density) as well as the light element abundances more 
precisely (see the lectures by Schramm and Walker). 

As a consequence of these observational advances, cosmology is becoming data- 
driven in an unprecedented way: theorists no longer have the luxury of untethered 
speculation, but must now confront their models with an impressive array of ob- 
servations. There is still debate about the reliability and interpretation of much 
of the data, but we are definitely entering the 'scientific' age of cosmology: those 
theories which are sufficiently worked out are becoming increasingly falsifiable and 
many will stand or fall in the coming years. This is a very healthy development 
for the field. It is safe to say that at present the big bang framework for the large- 
scale evolution of the universe remains healthy, but that we still lack a standard, 
tested model for the origin and evolution of structure within this framework (see 
the lectures by Scherrer on structure formation). 

In these introductory lectures, I cover only a very small portion of the many 
topics of recent interest in cosmology. The first chapters provide a general overview 
of the standard cosmology, the hot big bang model, focusing on its kinematics and 
dynamics, the observational evidence in its favor, and the current status of mea- 
surements of the global cosmological parameters. Subsequent chapters review the 
early universe, exploring aspects of the cosmic microwave background radiation and 
particle relics. The final section presents a brief introduction to the inflationary sce- 
nario for the very early universe and some of its observational implications. The 
inflationary scenario has not yet received the empirical backing to justify its inclu- 
sion in the 'standard' cosmological model; it is discussed here because it nevertheless 
provides a compelling theoretical framework for the very early universe, one which 
should be tested by observations in coming years. 

Readers wishing to delve more fully into these topics would do well to move 
on to two textbooks which cover them in substantial detail: Principles of Physical 
Cosmology, by P. J. E. Peebles 1 , and The Early Universe, by E. W. Kolb and M. S. 
Turner 2 . 

Before plunging in, I note that the appendix contains a brief discursion on 
notation and units. 
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2. The Standard Cosmology 



2.1. Homogeneity, Isotropy, and the Cosmological Principle 

The standard hot Big Bang model, based on the homogeneous and isotropic 
Friedmann-Robertson- Walker (FRW) spacetimes, is a remarkably successful oper- 
ating hypothesis describing the evolution of the Universe on the largest scales. It 
provides a framework for such observations as the Hubble law of recession of galax- 
ies, interpreted in terms of the expansion of the universe; the abundances of the light 
elements, in excellent agreement with the predictions of primordial nucleosynthesis; 
and the thermal spectrum and angular isotropy of the cosmic microwave background 
radiation (CMBR), as expected from a hot, dense early phase of expansion. 

While homogeneity and isotropy are, strictly speaking, assumptions of the model, 
they rest on a strong and growing foundation of observational support. The evi- 
dence for angular isotropy on large scales comes from the smallness of the CMBR 
large-angle anisotropy detected by COBE 3 and FIRS 4 (the quadrupole anisotropy 
detected in the first year of COBE data is (AT/T)/= 2 ~ 5x 1CT 6 ), from the isotropy 
of radiation backgrounds at other wavelengths, as well as from the isotropy of deep 
galaxy and radio source counts. Note that these different sources carry information 
about quite disparate scales and types of matter in the Universe. For example, 
the APM 5 and EDSGC 6 surveys measured the angular positions of more than one 
million galaxies over an area covering about 10 % of the southern sky out to an 
effective depth of roughly 600 h~ x Mpc. The galaxy surveys give us information 
about the local distribution of luminous matter in the universe. On the other hand, 
through the Sachs- Wolfe effect, the large- angle CMBR measurements directly probe 
the gravitational potential over length scales of order 6000 h~ x Mpc, and are thus 
sensitive to the mass distribution itself - the distinction is important to keep in 
mind, since the evidence for dark matter should lead us to be wary about identify- 
ing the distribution of light with that of mass. We will return to the CMBR below 
and first focus on the galaxy distribution. 

In a survey of galaxy angular positions to some limiting apparent brightness, 
the joint probability of finding two galaxies in elements of solid angle dfl\ and df^ 
separated by an angle 6 is given by 

dP = iV 2 dr2idr2 2 [l + w gg (6)\ , (2.1) 

where N is the mean surface density of galaxies in the survey and w gg (9), the 
galaxy two-point angular correlation function, measures the excess probability over 
random of finding a galaxy pair with this separation. If galaxies are distributed 
isotropically on large scales, we should find w gg (9) — » at large angles, and the 
correlation function should scale in a definite way with survey depth, both of which 
are observed. The APM data for w gg are shown in Fig. 1 (from 5 ). The correlation 
function is a power law, w gg (9) ~ #-°- 67 f or Q <; i° ? but breaks below this power law 
at 6 ~ 3°; for 6 <: 6°, \w gg {6)\ ^ 5x 10~ 4 and becomes lost in the noise. This behavior 



3 



0.1 - 



w(e) 



o.oi 



o.ooi 




e (deg) 



Figure 1: Galaxy angular correlation function w gg (9) measured in the APM survey, 
for galaxies in the apparent magnitude interval 17 < bj < 20. 



implies that the variance in the number of galaxies in a patch of fixed angular size 
becomes small for large patch size: in other words, averaged over sufficently large 
angular scales, we see roughly the same number of galaxies (brighter than a fixed 
limit) per steradian in different parts of the sky. 

Evidence for large-scale homogeneity comes in part from galaxy redshift surveys. 
However, compared to the angular photometric surveys, which currently give two- 
dimensional information for several million galaxies, currently complete redshift 
surveys yield three-dimensional information for typically several thousand galaxies 
in a more local neighborhood. In the redshift surveys, one can verify directly that 
the rms fluctuations in the spatial number density of galaxies become small when 
averaged over large enough scales. For example, in the full-sky surveys selected from 
infrared sources in the IRAS catalog (the 1.2 Jy survey of Fisher, etal. 7 and the 1- 
in-6 QDOT survey of Efstathiou, etal. 8 shown in Fig. 2), the rms fluctuation in the 
number of galaxies in cubical volumes of side L = 60/i — 1 Mpc is of order SN ga i /N ga i ~ 
0.2 and decreases with increasing cell volume. This approach to homogeneity is 
roughly consistent with that seen in the much deeper (but two-dimensional) angular 
surveys. [In fact, the approach to homogeneity as a function of increasing scale 
observed in Fig. 1 is somewhat more gradual than was expected in the popular 
cold dark matter model for galaxy formation-this is the famous problem of extra 
large-scale power, which we will comment upon later.] Larger structures such as 
superclusters, great attractors, great voids, and long filaments do exist, and have 
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Figure 2: An equal-area projection of the IRAS-QDOT galaxies (open circles) at 
distances in the range 20 — 500/i _1 Mpc. Solid squares are regions not surveyed, and 
the black band is an excluded region encompassing the Galactic disk (from ref.10). 

received considerable attention. Of particular note in this regard was the discovery 
of the Great Wall, extending roughly 170 X 60 X 5h~ 3 Mpc 3 , by Geller and Huchra 9 
in the Center for Astrophysics (CfA) survey extension. But in a statistical sense, 
the net fluctuations in galaxy number become small on the largest scales where they 
have been reliably counted in large-area redshift surveys. This is consistent with 
the visual impression from Fig. 2. 

This trend is confirmed by the visual appearance of large-scale structure in the 
ongoing Las Campanas redshift survey 11 : going considerably deeper than the CfA 
survey, they do not find evidence for coherent structures larger than the Great 
Wall. The Las Campanas survey uses a multifiber spectrograph to simultaneously 
measure many redshifts in the same field, an ingenious technological development 
which makes possible the extension of complete large- area redshift surveys to greater 
depths in a finite survey time. Advances in multi-fiber spectroscopy will be further 
exploited by the Sloan Digital Sky Survey, which aims to measure one million galaxy 
redshifts over a contiguous area of n sr in the northern sky. This survey will use 
a 600-fiber spectrograph to accumulate redshifts at an unprecedented rate. A con- 
current photometric survey will measure angular positions for roughly 50 million 
galaxies. 

From the CMBR and galaxy observations discussed above, we infer that on 
large scales the universe appears to be distributed in a statistically homogeneous 
and isotropic way around us. To step from this evidence to a cosmological model, 
we must add to it a further assumption, because we observe the universe from only 
a single vantage point and have direct information only about conditions on our 
past light cone. This assumption, often called the Copernican Principle, states that 
we do not occupy a special position in the Universe: the conditions we observe are 
typical of those seen by observers on distant galaxies. By itself, the Copernican 
Principle would allow a large range of permissible cosmological models, e.g., an 
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inhomogeneous fractal universe with self-similar structure on all scales, or homo- 
geneous models with anisotropic expansion in different directions. However, when 
combined with the observations of isotropy above, the Copernican principle becomes 
quite powerful, for it implies that the Universe should appear isotropic about every 
point. By a straightforward geometric argument, isotropy about every point (or 
more precisely, about every local fundamental observer comoving with the cosmic 
fluid) in turn implies that the Universe is spatially homogeneous, i.e., that there 
exist spacelike 3-surfaces of uniform energy-momentum which evolve according to a 
universal time t. This more general assumption of global homogeneity and isotropy 
is called the Cosmological Principle, a symmetry principle that is the foundation 
of the standard cosmology. 

At this point, one is naturally tempted to ask why the universe should have such 
strong symmetries of isotropy and homogeneity. The attempt to provide a dynam- 
ical answer to this question is one of the motivating forces behind the inflationary 
scenario, which we shall discuss later on. For now, we will adopt the cosmological 
principle as a working hypothesis and explore some of its consequences. 

2.2. Hubble's Law and Particle Kinematics 

The cosmological principle leads directly to Hubble's law of the expanding uni- 
verse. Consider a triangle formed by three fundamental observers (each of whom 
defines the local cosmic rest frame) at some initial time. If the universe remains 
homogeneous and isotropic over time, then the triangle at all later times must be 
similar to the original one: the length of each side scales up by the same factor a{t). 
Extending this argument to all other fundamental observers, we see that a(t) must 
be a universal scale factor, such that the distance £{t) between any two fundamental 
observers satisfies £(t) = £oa(t), where Iq is the initial separation. Then the relative 
speed of one observer with respect to the other is given by 

v(t) = d£/dt = h£ = (h/a)£(t) = H(t)£(t) . (2.2) 

Note that H is in general time-dependent, but for observers sufficiently nearby 
that the light travel time £jc is small compared to the time over which H changes 
appreciably, we can replace H by its present value, H(to) = Hq, where subscript 
denotes the present epoch. Relation (2.2) then reproduces Hubble's observation in 
1929 that the recession speed of a galaxy is proportional to its distance from us. 
The proportionality constant, Hq = (a/a)o, is now known as the Hubble parameter. 

A recent example of this is shown in Fig. 3, taken from ref. 12 : each circle indi- 
cates a cluster of galaxies, with distance estimated using the Tully-Fisher relation 
for spiral galaxies and recession velocity inferred from the redshift (see below); the 
horizontal error bars indicate the 1 — a spread in inferred distances for a number 
of galaxies per cluster. The slope of the straight line fit to the points is Hq ~ 80 
km/sec/Mpc; other methods of estimating Hq yield values in the range 40 - 100 
km/sec/Mpc, and it is conventional to write this as Hq = lOOh km/sec/Mpc, with 
0.4 < h < 1. (The closed circles around 40 Mpc indicate clusters thought to be 
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Figure 3: Hubble's law: galaxy recession speed is shown against distance inferred 
from the Tully-Fisher relation. 



falling into the Great Attractor region, an inhomogeneity that locally perturbs the 
Hubble flow.) 

The cosmological principle, with eqn.(2.2), tells us about the kinematics of par- 
ticle motions in the universe. First consider non-relativistic particles. A massive 
free particle passes a fundamental observer FOi at time t with relative speed v p <C c: 
since the FO is at rest with respect to the local fluid, v p is the peculiar velocity of 
the particle with respect to the local rest frame. After a time interval dt, the par- 
ticle has moved a distance di = v p dt and overtakes a second fundamental observer 
F0 2 , who has a speed dv = H(t)dl = Hv p dt relative to FOi. At that time, F0 2 
measures the particle's peculiar velocity to be v p (t + dt) = v p (t) — dv. Thus, the 
peculiar velocity satisfies the equation of motion 



dv p 
~dt 



dv 
"dt 



The solution is 



V p (X 



(2.3) 



(2.4) 



that is, in a homogeneous and isotropic universe, peculiar velocities decrease with 
the expansion. Now consider a dilute non-relativistic gas in thermal equilibrium. 
By equipartition, the mean kinetic energy is m(v 2 )/2 = (3/2)/c£T fl , where T g is 
the gas temperature; by eqn.(2.4) the gas temperature satisfies T g ~ v 2 ~ a~ 2 . 
This applies, for example, to baryons well after they decouple from the photon 
background (before decoupling, the baryon temperature is tied by scattering to the 
photon temperature, which scales differently-see below). 
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The same reasoning can be applied to study the evolution of relativistic particles. 
Consider a photon with initial frequency u(t) as measured by FOi. After the interval 
dt, it has travelled a distance di = cdt and passes FO2, who has a speed v = Hd£ = 
Hcdt with respect to FOi. The photon frequency v(t + dt) observed by FO2 is given 
by the first-order Doppler shift, 

dv = u(t + dt) - u(t) = -v- = -uHdt . (2.5) 

c 



The frequency thus satisfies 



with solution 



dv a , . 

-r = -v- , 2.6 

dt a 

u(t)cc^-r . (2.7) 
a{t) 

The photon frequency is redshifted with the expansion, and its wavelength is stretched 
with the scale factor, \(t) ~ v~ x ~ a. For a photon gas in thermal equilibrium, 
the temperature thus scales inversely with the scale factor, T r ~ (E r ) ~ v ~ a -1 . 
Generalizing the argument, one finds that the De Broglie wavelengths of all free 
particles follow the scale factor in this way. This leads us to define the redshift z: 

l + zM = ^ = a JM (2.8) 

is the ratio of the photon wavelength at emission (t e ) to its wavelength observed at 
the present (to). This expression holds more generally than this derivation suggests: 
it is valid even when the period v~ x is comparable to or longer than t (although 
its operational meaning is then less transparent). The redshift z thus plays several 
roles: through the Doppler shift, it is a measure of recession velocity, v ~ cz; through 
the Hubble law, it is a measure of distance, v = cz = Hq£; and through eqn.(2.8), 
the redshift z(t) can be used to characterize a cosmological epoch t. The most 
distant objects directly observed are luminous quasars, which have been seen out to 
redshifts approaching z = 5. By comparison, the photons in the cosmic microwave 
background radiation were probably emitted (or, more precisely, last scattered) at 
a redshift z ~ 1000. 

2.3. The Metric: Friedmann-Robertson- Walker Models 

In the context of general relativity, the cosmological principle severely restricts 
the form of the spatial geometry: a uniform stress-energy momentum tensor implies 
that the constant-time 3-surfaces have uniform spatial curvature. The metric on 
such surfaces has the form 



, 9 9 , 

dsn = a it) ■ — : 

3 W 1 - jfcr 



+ r 2 (de 2 + sin 2 9d<f) 



(2.9) 
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where a(t) is the global scale factor which describes the overall expansion or con- 
traction, r,9,(f) are the (fixed) comoving coordinates carried by the fundamental 
observers, and k = 0, 1, —1 is the sign of the spatial curvature. The case k = cor- 
responds to flat, Euclidean 3-space (i? 3 ) written in spherical coordinates, k = 1 cor- 
responds to the geometry of the three-sphere (S 3 ), and k = — 1 to the 3-hyperboloid 
(H 3 ), the three-dimensional analogue of a hyperbolic saddle. Thus, models with 
k < are spatially infinite (open), while those with k = 1 are spatially finite 
(closed). (Note that the apparent singularity at r = ±1 for k = 1 is the usual 
coordinate singularity at the poles of a sphere.) By homogeneity and isotropy, the 
full spacetime metric then takes the Friedmann-Robertson- Walker (FRW) form 

ds 2 = g™dx^dx" = dt 2 - ds\ , (2.10) 

where t is proper time measured on the clocks carried by the fundamental comoving 
observers. 

We can now be somewhat more precise about the assumption of homogeneity and 
isotropy: clearly this is meant to be taken in some average sense over large scales. It 
is useful to think of this perturbatively: the actual spacetime metric can be written 
as the FRW metric plus a perturbation: g^ v = g^ w + /i M „(x,t). In an appropriate 
gauge, /ioo satisfies an equation of motion analogous to the Newtonian potential 
6. Through the Sachs- Wolfe effect, the COBE observations roughly indicate that 
ST/T ~ <fi ~ /iqo ~ 10~ 5 <C gQQ RW = 1: the departure from the FRW metric on large 
scales is very small (at least in the chosen gauge). This argument extends to smaller 
scales as well: although galaxies represent highly non-linear condensations of the 
mass density, their associated gravitational potential is small, <f> ^ 10~ 4 . Thus at 
least when averaged over scales larger than individual galaxies, the perturbation to 
the FRW spacetime metric associated with inhomogeneities is very small, and the 
homogeneity assumption is an excellent first approximation. 

The FRW models are characterized by the global scale factor a(t), whose dynam- 
ics is determined by the matter content of the universe through Einstein's equations, 

~ \g^R - A9n„ = 8nGT^ . (2.11) 

Treating the matter as a perfect fluid locally at rest, and using the FRW metric, 
these become the Friedmann equations 

H 2 = - =— P--J + - 2.12 

\aj o cr 6 

and 

a = _4*G(p + 3p) A 

a 3 3 V ' 

Here p = Tqq is the mean energy density of matter, p = Ta is its pressure, and A is the 
cosmological constant, i.e., the effective contribution to the energy-momentum from 
the vacuum state. Although it follows from the Einstein equations, it is also useful 
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to separately write down the energy-momentum conservation equation, V v T jlv = 0, 
for the fluid in the FRW universe, 

^ + 3H(p + p) = . (2.14) 

Observations suggest that the fluid energy density of the universe is currently 
dominated by non-relativistic matter (m), while the early universe was dominated 
by ultrarelativistic particles, or radiation (r). Non-relativistic matter is also some- 
times called pressureless dust, because its fluid pressure is dynamically negligible, 

Pm ~ Pm(v 2 ) < p m ■ ( 2 - 15 ) 

In this case, from eqn.(2.14), the energy density scales as 

Pm ~ a" 3 ■ (2.16) 

This is consistent with conservation of particle number in a fixed comoving volume 
(n m ~ a -3 ), since the particle energy is dominated by its rest-mass, p m = mn m . For 
radiation, the equation of state is 

Pr = Pr/3 , (2-17) 

and eqn.(2.14) then yields 

p r ~ a~ 4 . (2.18) 

The faster fall-off of p r with the scale factor is consistent with the Doppler redshift of 
the frequency of a photon emitted by a receding observer in the expanding universe, 
eqn.(2.7). For massless particles with average energy (E r ), this implies (E r ) ~ a -1 
and thus p r ~ (E r )n r ~ a~ 4 . 

Two features of the solutions to the Friedmann equation (2.12) with A = 
are worth noting. In this case, there is a one-to-one correspondence between the 
spatial geometry and the fate of the universe: open models (A; < 0) expand forever, 
while closed models (k > 0) eventually recollapse, because the energy density for 
matter and radiation fall off faster than a~ 2 . Conversely, in the early universe, 
a <C «o, the matter and radiation terms dominate over the spatial curvature, and 
the dynamics of the model is well approximated by setting k = 0. In this limit, for 
a matter-dominated universe, eqns.(2.12) and (2.16) give the solution 

a(<)~< 2/3 ; Pra = ^r 2 (MD,& = 0) (2.19) 
while for radiation-domination, from (2.12) and (2.18) 

a(i)~< 1/2 ; Pr = ^^ (RD,A; = 0) . (2.20) 
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This clearly generalizes for a fluid wifh equation of state pj = wpf. 

-3(l+«j) _^ r /.L\ + 2l2(l+w) . „ \ /n ni\ 

Pf~ a ( >^a(t)~t" , Pf=^ Gt2{l + wf ■ (2-21) 

In the next few sections we focus on the cosmologically recent matter-dominated 
era, and subsequently turn to the radiation-dominated early universe. 

3. Cosmological Parameters: Hq and to 

The principal observable cosmological parameters of the FRW models are the 
Hubble parameter, Hq = (a/a)o, the age of the Universe, tQ, the present (non- 
relativistic) mass density relative to the 'critical' density of the spatially flat, Einstein- 
de Sitter (A; = A = 0) model, 

O = Po/pcrit = 87rGp /3i?o , (3.1) 

the deceleration parameter, qq = — (aa/a 2 )o, which measures the rate at which the 
gravitational attraction of the matter is slowing down the expansion, and the con- 
tribution of the cosmological constant to the present expansion rate, Ao = A/3Hq. 
From the Friedmann equations, these parameters are related by 

1 = + Ao - -^72 ( 3 -2) 

For vanishing cosmological constant, fig — 1 determines the sign of the spatial cur- 
vature: Qq = 1 for the spatially flat model (k = 0), and it is less than one for open 
models. 

The age of the universe is related to the other parameters through an expression 
of the form i?oto = 1.02/z(to/10 10 yr) = /(Jlo,Ao), where / is a function of order 
unity. Thus the Hubble time Hq 1 = 9.8 X 10 9 /i _1 yr sets the timescale for the age of 
the universe, while the Hubble length cHq 1 = 3000/i _1 Mpc sets the lengthscale for 
the present observable universe. For a matter-dominated universe with A = 0, / 
falls monotonically with increasing Qq, and two useful limits are /(0,0) = 1 (a ~ t 
for Q,q = 0) and /(l, 0) = 2/3 (since a ~ t 2 ^ 3 for O = 1). More generally, over the 
range < Qq < 1, k < 0, O — 3Aq/7 < 1, an excellent approximation is 13 



2 sinh- 1 ( v /(l- OO/O,) 

m ° * 3 — TT^m — (3 - 4) 

where 

O = O-0.3(O + A ) + 0.3 . (3.5) 
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Figure 4: Evolution of the cosmic scale factor (scaled to its present value). 

Some examples of the evolution of the scale factor a(t) are shown in Fig. 4, 
including three cases with A = 0: empty (Qq = 0), flat (Qq = 1), and closed 
(Qq = 5), and one example of a flat (k = 0) model dominated by a cosmological 
constant (Aq = 0.95 = 1 — f^o)- This figure displays the decrease in the expansion age 
HqIq from unity as Qq is increased from zero. It also demonstrates that a significant 
cosmological constant term can make HqIq large. The solution for the spatially flat 
(k = 0) model with non-zero A has the form 



where Aq + S1q = 1. At late times, it approaches the exponential de Sitter solution, 



Much early effort was spent trying to measure or constrain the parameters q$ 
and Ao through the classical 'cosmological tests', such as the Hubble diagram, an- 
gular size as a function of redshift, and galaxy counts as a function of redshift and 
apparent brightness. For example, to construct the Hubble diagram, one measures 
the apparent brightness of a well-defined sample of objects (say, the brightest galax- 
ies in clusters) as a function of the object's redshift; for galaxies of fixed intrinsic 
luminosity, the scaling of apparent magnitude with redshift is a function of the 





a 



(3.7) 
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cosmological parameters. Unfortunately, galaxies at large distances, where the dis- 
tinction between model parameters becomes observable, are seen when they were 
much younger than their nearby counterparts, so a model for galaxy luminosity 
evolution must be used to interpret the results. Significant progress has been made 
in understanding galaxy evolution, and there is hope that the effects of evolution 
and cosmology might be disentangled in coming years, but these tests currently do 
not place stringent universally accepted constraints on the cosmological parameters. 
Recently, it has been pointed out 14 that the probability that a quasar at a given 
redshift is gravitationally lensed by a foreground galaxy is a sensitive test for the 
cosmological constant: in models with Aq > 0, one generally expects a higher lens- 
ing probability Based on surveys for lensed quasars, the bound Ao 0.8 has been 
inferred in the case of a spatially flat (k = 0) universe 15 . For more on the standard 
cosmological tests, the reader is referred to Peebles 1 . 

3.1. Hq and the Distance Scale 

The Hubble parameter relates the observed recession velocity v r or redshift z of 
a galaxy to its distance d: for v r . <C c, the recession velocity is v r . = cz = H^d + v p , 
where v p is the peculiar radial velocity of the galaxy with respect to the Hubble 
flow, usually assumed to arise from gravitational clustering. Galaxy redshifts can be 
measured quite accurately, so all the difficulty in determining Hq resides in finding 
reliable distance indicators for extragalactic objects at distances large enough that 
the Hubble term dominates over the peculiar motion. Observed peculiar velocities 
are typically of order 300 km/sec, so that distance measurements beyond 40 Mpc or 
more (recession velocities above 4000 km/sec) are required for reasonable accuracy. 

A wide variety of techniques has been used to establish an extragalactic dis- 
tance scale 16 , and this is reflected in the spread of results for Hq, roughly 40 — 100 
km/sec/Mpc. Distance estimates made using methods such as the Tully-Fisher 
relation between 21-cm rotation speed and infrared luminosity for spiral galaxies, 
calibrated by observations of Cepheid variable stars in several nearby galaxies, have 
yielded high values for the expansion rate, roughly Hq = 80 ± 10 km/sec/Mpc. Two 
newer methods, planetary nebula luminosity functions 17 and galaxy surface bright- 
ness fluctuations 18 yield values for Hq in this range as well, and are being further 
developed. On the other hand, methods using Type la supernovae as standard can- 
dles have yielded low values, Hq ~ 50 ± 10 km/sec/Mpc. SNe la are thought to be 
the explosions of white dwarfs which accrete matter from binary companions until 
they reach the Chandrasekhar mass, and there is some evidence that they form a 
homogeneous class; they also have the advantage that they can be observed to great 
distances. In the future, Hubble Space Telescope observations of Cepheids in other 
nearby galaxies (and perhaps as far as the Virgo cluster) which are hosts to SNe la 
or which can be used as Tully-Fisher calibrators should help improve the situation. 
The recent discovery of a probable Type la supernova at z = 0.46 19 also raises hopes 
that a sample of SNe la at z ~ 0.5 could significantly constrain qq, provided the 
dispersion in SNe la luminosities is sufficiently narrow. 
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There are also a variety of methods being employed to measure the distances of 
extragalactic objects directly, bypassing the extragalactic distance ladder built up 
from Cepheids. Using the expanding photosphere method, Schmidt, Kirshner, and 
Eastman 20 have determined the distances to 10 type II supernovae at large distances, 
and find good agreement with the Tully-Fisher distances for these galaxies. Other 
'direct' methods which hold future promise include measurement of the Sunyaev- 
Zel'dovich effect, due to the Compton upscattering of CMBR photons by hot gas 
in rich clusters 21 , and the differential time delay between images in gravitationally 
lensed quasars 22 . 

3.2. The Age of the Universe, t 

Three methods have traditionally been used to infer the age of the Universe, 
to- Nuclear cosmochronology is based on radioactive dating of r— process elements, 
that is, heavy elements formed by rapid neutron capture, most probably in super- 
novae. The element ratios Re/Os and Ur/Th generally indicate to = 10 — 20 Gyr 23 , 
with a large uncertainty due to the unknown element formation history (e.g., the 
star formation rate over time). A second method involves the cooling of white 
dwarfs: when low-mass stars exhaust their nuclear fuel, they become degenerate 
white dwarfs, gradually cooling and becoming fainter. The number of white dwarfs 
as a function of luminosity drops dramatically for L w d < 3 x 10~ 5 L Q , suggesting 
that there has not been sufficient time for them to fade below this value. Coupled 
with models of white dwarf cooling, this implies that the age of the galactic disk is 
about t ^ 10 ± 2 Gyr 24 . 

The most extensively studied technique for constraining to i s the- determina- 
tion of the ages of the oldest globular clusters in the galaxy. When stars finish 
burning hydrogen, they turn off the main sequence, characteristically reddening 
and brightening. By observing the color-magnitude (color vs. apparent bright- 
ness) diagram for a cluster, one can determine the apparent brightness of stars 
in the cluster that are now leaving the main sequence. Knowing the distance to 
the cluster then gives the absolute luminosity of stars at the turn-off. On the 
other hand, stellar evolution theory relates stellar luminosity to the time a star 
spends on the main sequence: crudely, a star turns off the main sequence when 
a fixed fraction / ~ 0.15 of its initial hydrogen mass Mh = XhM has been con- 
verted into helium (here Xh ~ 0.75 is the initial H mass fraction of the star). 
The conversion of H to He releases 6.4 x 10 18 erg for every gm of H burned. 
Thus, over the course of the main sequence, the total energy radiated by a star 
is E oc fXjjM ~ 10 52 /X# (M/Mq) erg, and its evolutionary lifetime on the main 
sequence is r ms ~ EjL ~ 10 11 fXjj(M /M Q ) /(L/ ' L & ) yr. Now, stellar models indi- 
cate that the relationship between mass and main sequence luminosity is roughly 
L ~ M 4 , yielding a main-sequence lifetime of T ms ~ 1.2 X 10 10 (L/ J L Q ) _3 ' /4 years. 
Thus, the measurement of the main sequence turn-off luminosity gives an estimate 
of the age of the cluster. The turn-off luminosity in the oldest globular clusters in the 
galaxy is slightly below that of the sun, yielding the age estimate t gc = (13 — 15) ±3 
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Gyr . This argument is a gross oversimplification of the complex process of stellar 
modelling and cluster isochrone fitting, but it gives a heuristic understanding of the 
resulting age estimate. The largest source of error is apparently the uncertainty 
in the distances to the globular clusters. It is hoped that observations with the 
corrected Hubble Space Telescope mirror, successfully repaired in late 1993, could 
reduce the uncertainty in t gc to as little 10%. There may also be residual systematic 
model uncertainties associated with the initial hydrogen (or helium) fraction Xh, 
opacities, metallicity, etc. 

We now turn to the final cosmological parameter of interest, the density of the 
universe. 

4. Cosmological Parameters: Qo and Dark Matter 

It is convenient to parameterize the mass density of the universe in terms of 
the mass-to-light ratio, say in the Vr band, T = (M/L)/(Mq/L q ). Dividing the 
present critical density, 

p c = 3H%/8ttG = 1.88/z 2 x l(T 29 gm cnT 3 = 2.8 x 10 n h 2 M Q Mpc~ 3 (4.1) 

by the observed mean luminosity density jy T ~ 2.4 x 10 8 hL Q Mpc~ 3 , the critical 
mass-to-light ratio for the Qq = 1 universe is T c ~ 1200/i, and the cosmic density 
parameter can be expressed as H = 8 x 1CT 4 /i _1 T. The mass-to-light ratio in 
the solar neighborhood is approximately T = 5, while the central cores of ellipti- 
cal galaxies yield T ~ 12h, so the density of luminous matter (that is, of matter 
associated with typical stellar populations) is inferred to be fli um ~ 0.007. How- 
ever, it is well known that the luminous parts of galaxies are not the whole story: 
there is strong evidence from flat spiral galaxy rotation curves, from a variety of 
observations of galaxy clusters, and from large-scale peculiar motions that there is 
a substantial amount of dark matter in the universe. 

4.1. Dark matter in galaxies 

An example of this phenomenon for galaxies is shown in Fig. 5, which shows a 
model for the rotation speed of the galaxy as a function of radius. The three lower 
curves show the contributions from the luminous components (two bulge compo- 
nents and the disk 26 ). Clearly, they cannot account for the observed rotation speed 
beyond a few kpc from the galactic center. To match the observations (taken from 
the compilation in ref. 27 ), one assumes in addition a quasi-spherical distribution of 
dark matter, called the halo, with a density profile of the form 



In Fig. 5, the adopted halo core radius is a = 3.5 kpc. On large scales, the halo mass 
scales as Mh(r) ~ r, yielding a flat rotation curve, w 2 = GM(r)/r — >• 220 km/sec 
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Figure 5: Rotation curve for our galaxy: an average over a number of different 
observations (solid points) is shown with a model that includes luminous and dark 
matter. 

for the model shown in Fig. 5. Clearly, eqn.(4.2) cannot extend to arbitrarily large 
radii, because the enclosed mass would diverge; the halo is presumably truncated by 
tidal or other effects at large r. The observation of high proper motion stars in the 
solar neighborhood (presumably bound to the Galaxy) implies that the local value 
of the galactic escape velocity exceeds 450-500 km/sec, and this can be used to 
place a lower bound on the halo truncation radius. If the halo is sharply truncated 
at r t , the escape speed v e at r < r t satisfies 28 



Using v e > 500 km/sec at the solar radius r = Rq = 8.5 kpc and v c = 220 yields 
r t ^ 40 kpc and a total galaxy mass M ga \ ^ 5 X 10 11 M Q . Comparing with the 
luminosity of the disk plus the bulge, L ga j = 1.4 x 10 10 L Q in the V band, this 
implies that the total mass-to-light ratio for the Milky Way is at least Tmw ^ 35. 
This is consistent with the requirement that distant globular clusters and satellite 
galaxies are bound to the Galaxy, as well as with mass-to-light ratios inferred from 
more extended flat rotation curves in other spiral galaxies 27 . If these systems are 
typical of the universe, we infer Qq <: 0.02/j _1 for the matter associated with galaxy 
halos. 

It is interesting to compare these values with the baryon density inferred from 
primordial nucleosynthesis (see the lectures by Schramm and Walker), which 29 30 
has been restricted to the range 0.010 < f^/z 2 < 0.015. Comparison with fli um 
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above suggests that some or most of the baryons are dark or in underluminous 
populations. Furthermore, comparison with the escape speed and rotation curve 
bound Qq <: 0.02/j _1 shows that baryons could constitute some or all of the dark 
matter in galaxy halos, depending on their extent. 

One possibility for dark baryons in halos would be degenerate brown dwarfs, 
substellar (M < O.O8M ) objects which did not reach sufficiently high tempera- 
ture to burn hydrogen. Currently several independent groups are searching for halo 
dwarfs (which have been dubbed MACHOS, for massive compact halo objects); the 
signature is a microlensing event, in which a background star, say in the LMC or the 
bulge of the Milky Way, symmetrically brightens and fades as a MACHO passes 
near its line of sight. 31 Although they are distinguishable from ordinary variable 
and flare stars by the time symmetry and achromaticity of the light curve, such 
microlensing events would be intrinsically rare-the probability that a given back- 
ground star is lensed is r ~ (v c /c) 2 ~ 5 x 10~ 7 , so a large number of stars must 
be accurately monitored. The American- Australian MACHO project currently ob- 
tains CCD photometry for several million stars per night, and has discovered many 
periodic variable stars. The French EROS and Polish OGLE groups are also mak- 
ing progress. In September and October 1993, all three groups reported the first 
evidence for microlensing events by objects in the Galaxy 32,33 ' 34 . The most prob- 
able masses for the lensing objects are roughly in the range 0.01 — O.5M . As 
data is accumulated over the next few years, we will be able to infer information 
about the contribution of MACHOs to the Milky Way halo. It is also worth noting 
that microlensing in external galaxies has been sought previously by monitoring 
the brightness of gravitationally lensed QSO images over time. In particular, a mi- 
crolensing event in the lensed quasar 2237+0305 has been established, 35 but in this 
case it is difficult to establish the MACHO mass, since the microlensing probability 
is closer to unity. 

4.2. Dark matter in clusters 

Moving to larger scales, for clusters of galaxies the traditional dynamical method 
of estimating cluster masses and mass-to-light ratios, first used by Zwicky who 
discovered the 'missing mass' problem in the 1930's, has been to apply the virial 
theorem to measured cluster velocity dispersions, M tot oc {v 2 )/{Rjj), where {v 2 ) 
is the velocity dispersion of the galaxies in a cluster and R{j is the separation 
between them. This method assumes that galaxies trace the cluster mass and 
that the galaxy velocity distribution is isotropic, both of which may not be good 
approximations. Independent information on the dark matter distribution in the 
inner cores of clusters comes from the giant luminous arcs and arclets, high redshift 
galaxies gravitationally lensed by foreground clusters 36 37 . These arcs are formed 
when a galaxy is nearly imaged into an Einstein ring. Measurements of the cluster 
and background galaxy redshifts yield an estimate of the cluster mass within the 
impact parameter of the lens; for most of the cases studied, these estimates are in 
reasonable agreement with the mass-to-light ratios inferred from the virial theorem, 
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but it should be noted that the arc observations probe only the inner few hundred 
kpc of the clusters. (Work is presently being done to extend this idea to map clusters 
on larger scales, by looking for more minute statistical distortions in the alignments 
of background galaxies.) Typical inferred cluster values are T ~ 100 — 250/i, which 
would imply Q$ ~ 0.1 — 0.2. Since clusters are rare objects, occupying a very small 
fraction of the universe, one expects this to be a lower bound, in which case some 
form of non-baryonic dark matter is probably required, given the limits on Qb above. 

X-ray observations, most recently by the ROSAT satellite, 38 have begun to map 
the density and temperature profiles of the hot gas which permeates many clusters. 
Since the gas is in hydrostatic equilibrium, it can be used to trace the cluster mass 
distribution (including dark matter) directly, 



where /i is the mean molecular weight of the gas, and n(r) and T(r) are the gas 
density and temperature. Cluster masses inferred from X-ray observations are gen- 
erally comparable to but ~30% less than virial estimates. Consistent with this, 
the X-ray measurements indicate that clusters are surprisingly baryon-rich, in the 
sense that the gas constitutes typically 5 — 10/i~ 3 ' 2 % of the inferred binding mass 
within approximately l/i _1 Mpc of the center of a rich cluster like Coma. (More- 
over, the dark matter tends to be more centrally concentrated than the gas out to 
this scale, suggesting that the gas fraction of the whole cluster is at least as large 
as the value above. 40 ) If this ratio is representative of the baryon mass fraction of 
the universe, then the nucleosynthesis bound on Qb would imply 41 the upper limit 
f^o ^ 0.15/i _1//2 . (Including the baryon stellar component in cluster galaxies only 
strengthens this limit.) This has been taken as evidence against the universe having 
closure density (Qq = 1) and would require advocates of inflation (which implies 
k = 0) to fall back on a cosmological constant or some other smoothly distributed 
component. The other possibility would be to loosen the nucleosynthesis bounds 
on Qb through some non-standard scenario such as inhomogeneous nucleosynthesis, 
but the upper bound on Qb is n °l raised sufficiently in this model to get around 
the argument above. One should perhaps be cautious about inferring the universal 
baryon fraction from a rich cluster like Coma. In particular, clusters with higher 
X-ray temperatures (and therefore larger total masses) may have a larger fraction 
of their total mass in gas. 40 This trend, coupled with the steeply falling distribution 
function of cluster temperatures, dn c /dT ~ T -5 for 3 < kT < 10 keV, suggests that 
the mean gas (and baryon) fraction may be substantially below the value for Coma, 
since the mean is dominated by the more numerous cooler clusters {kTc om a — 7 
keV). This would raise the derived upper bound on Qq closer to unity and suggests 
that massive clusters like Coma may not be representative of the baryon fraction 
of the Universe. On the other hand, N-body simulations of cold dark matter with 
baryons 41 indicate that the baryon fraction in a Coma-size cluster should be repre- 
sentative of the mean, so the present situation is somewhat confusing. 




(4.4) 
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4.3. Large-scale flows 

Moving to still larger scales, the deviations from the Hubble flow have been 
used to infer the cosmic density over scales of order 50/i _1 Mpc. The basic idea is 
to compare samples of the density perturbation field and the peculiar velocity field 
covering the same volume; 42 assuming they arise gravitationally, the proportionality 
between them depends on the rate of growth of the density fluctuations, which in 
turn depends on Qq. On very large scales, as discussed above, the rms density 
fluctuations are small, so linear perturbation theory away from the FRW spacetime 
is a reasonable first approximation (see the lectures by Scherrer). In this case one 
finds the relation 

V-v p = -H n° 6 8 , (4.5) 

where the density field 5(x) = (p(x) — p)/p, and the perturbation growth rate 
enters through d\n8/d\na ~ Q 6 . If one expresses distances d in terms of their 
equivalent Hubble velocities, v = H^d, then Hq drops out of Eqn.(4.5), so the 
uncertainty in the Hubble parameter does not undermine this method. A number 
of different approaches have been used to extract J7q in this way, and they have 
all given consistently high answers, Qq ^ 0.5. The topographic correlation between 
the observed galaxy field and the density field inferred from galaxy peculiar velocity 
surveys suggests that galaxies do broadly trace the mass distribution on large scales, 
in the sense that more galaxies are found in regions of high mass density (that is, 
high V • v p ), but the galaxy distribution may be 'biased' with respect to the mass. 
In the simplest linear bias model, the smoothed galaxy and mass density fields 
are assumed to be proportional, 8 ga i(x) = b ga iS(x), where b ga i is the bias factor, 
taken to be constant for a given class of galaxies. Since 8 ga i is what is observed, 
the comparisons based on Eqn.(4.5) actually constrain the combination il[j' 6 /6 fla /, 
where the bias factor refers, e.g., to galaxies selected from the IRAS catalog. Recent 
determinations have found fl^ 6 /b ga i ~ 1. For a bias factor of order unity, this is 
consistent with Qq = 1, which is pleasing to theorists and also buttresses the case 
for non-baryonic dark matter. It is well to keep in mind, however, that biasing 
is presumably a complex process associated with the non-linear stages of galaxy 
formation, so that the proportionality of the galaxy and density fields may be non- 
linear and/or scale-dependent. In addition, even with substantial smoothing of 
the density field, the perturbation amplitude in many regions is not small, and 
corrections to the linear theory must be taken into account. 

A summary of typical estimates of the density parameter on different scales is 
given in Table 1. 

4.4. Non-baryonic Dark Matter 

If CIq > Qb, as the dynamical observations on scales larger than clusters sug- 
gest, some form of non-baryonic matter must be invoked to make up the balance. 
Moreover, since this matter must be dark (Qq ^> Clj um ), it is natural to consider 
weakly interacting particles as candidates. It is convenient to distinguish two broad 
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Table 1: Estimates of the cosmic density parameter Qq. 



System 


T 


^0 


Luminous matter in galaxies 


5 - 1 2h 


0.007 


Galaxy halos 


35 


0.02/z- 1 


Clusters 


250h 


0.2 


Large-scale flows 




0.5 -f 



classes of non-baryonic dark matter, hot and cold, on the basis of their clustering 
properties. The prototypical hot dark matter candidate is a light neutrino with 
mass m v ~ 20 eV (see Sec. 6.3). Since they are relativistic {m v ^ T) until rela- 
tively recent epochs, light neutrinos would free-stream out of and damp out density 
perturbations up to the scale of galaxy clusters. Galaxies would form after clusters 
via fragmentation ('top down'). Due to phase-space constraints 43 , light neutrinos 
would not cluster significantly on the scale of galaxies: baryons would constitute 
the predominant dark matter in galaxy halos, while neutrinos would dominate in 
clusters. Cold dark matter, on the other hand, is defined to have negligible free- 
streaming length-it clusters on all scales. In cold dark matter models, structure 
generally forms hierarchically, with smaller clumps merging to form larger ones 
('bottom up'). For this reason, many (but not all) theorists since the early '80's 
have tended to prefer cold over hot dark matter, but it is worth noting the recently 
surging popularity of a mix 'n match scenario: a combination of 70% cold and 30% 
hot dark matter (say, with m„ ~ 7 eV) produces a favorable spectrum of large-scale 
density perturbations in the context of inflation, with some apparent advantages 
over pure cold dark matter. 

The theoretically favorite candidates for cold dark matter are weakly interacting 
massive particles (WIMPs). with masses generally in the range 20 — 150 GeV, and the 
axion, an ultra-light pseudoscalar with a mass of order 10~ 5 eV. The most attractive 
WIMP candidate is the neutralino, the lightest supersymmetric fermionic partner of 
the standard model bosons; its weak annihilation rate in the early universe naturally 
leaves it with an abundance comparable to the present critical density. The axion 
is the pseudo-Nambu-Goldstone boson associated with spontaneous breakdown of 
a global £7(1) symmetry (the Peccei-Quinn symmetry) introduced to explain why 
the strong interactions conserve CP. The global symmetry is spontaneously broken 
at some large mass scale fpQ, through the vacuum expectation value of a complex 
scalar field, ($) = fpg exp(ia/ /pq) / V%- At energies below the scale fpQ, the only 
relevant degree of freedom is the massless axion field a, the angular mode around 
the bottom of the $ potential. At a much lower energy scale, Aqcd ~ 100 MeV, the 
symmetry is explicitly broken when QCD becomes strong, and the axion obtains a 
periodic potential of height ~ Aqcd- ^ n 'invisible' axion models with Peccei-Quinn 
scale Jpq ~ 10 12 GeV, the resulting axion mass is m a ~ ^qcdI fpQ ~ 10~ 5 eV 
and fl a ~ 1. Although light, invisible axions interact so weakly, with cross-section 
a ~ 1//pq, that they were never in thermal equilibrium: they form as a cold Bose 
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condensate. 

Accelerator searches for neutrino mixing and beta-decay experiments on neu- 
trino mass should provide useful constraints on the possibility of neutrino dark 
matter. Active experimental efforts are also underway to detect both WIMPs and 
axions. Direct WIMP detection looks for the signals produced when a halo WIMP 
collides with a nucleus in a kg-size cryogenic crystal, depositing of order 10 keV 
in ionization and phonons (detector schemes based on scintillation and excitations 
of superfluids and superconductors are also being developed). Indirect WIMP de- 
tectors search for high energy neutrinos produced when WIMPs annihilate in the 
Sun and the Earth; large underground or underwater detectors currently in place 
or under development with sensitivity to WIMP annihilations include Kamiokande, 
MACRO, AMANDA, and DUMAND. Accelerator searches for supersymmetry also 
will constrain the neutralino parameter space. A large-scale axion search based at 
Livermore, expected to come on-line in 1994, will search for resonant conversion of 
halo axions to microwave photons in a cavity with a strong magnetic field. A scaled- 
up version of an idea originally proposed by Sikivie, this detector should approach 
the cosmologically interesting region of axion couplings for the first time 44 . 

Finally, given the paucity of direct evidence for dark matter, it is probably 
healthy to keep an open mind to alternatives. While it is natural to ascribe flat 
galaxy rotation curves and large cluster velocity dispersions to unseen matter, 
Milgrom 45 and others have argued that they may instead signal a breakdown of 
Newton's law of inertia at very low acceleration. The extent to which Milgrom's 
modified Newtonian dynamics (MOND) accounts for all the phenomena normally 
imputed to dark matter is controversial, and a full theory with which one could 
explore cosmology has been lacking. Nevertheless, at a minimum it provides a use- 
ful challenge to the accepted dogma. On a relatively more conservative side, the 
possibility that the dark matter interacts by other long-range forces in addition to 
gravity has recently been explored. 63 Such interactions are significantly constrained 
by galaxy and cluster observations, but could nevertheless have interesting impli- 
cations for structure formation and biasing. 

5. Cosmological Parameters: Taking Stock 

With this brief survey in hand, it is useful to pause and place these numbers for 
the cosmological parameters in theoretical context. If an extended period of inflation 
took place in the early universe (see below), then the spatial geometry should now be 
observationally indistinguishable from k = 0. If the cosmological constant vanishes, 
from Eqn.(3.2) spatial flatness implies fig = 1 (with the concomitant requirement 
of non-baryonic dark matter), and thus to = 2/3i?o = 6.5/i _1 x 10 9 yr. This is 
uncomfortably low compared to globular cluster ages unless h ^ 0.5 (to ^ 13 Gyr) 
and definitely problematic unless h < 0.65 (to > 10 Gyr), still on the low side of the 
Ho observations. However, a non-vanishing Ao is certainly allowed at some level by 
observations and has sporadically come into vogue, most recently to alleviate both 
this age problem and the large-scale power problem for inflationary cosmology. For 
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example, in a flat model with Qq = 0.25, Ao = 0.75, Eqn.(3.4) gives to ~ l/i?o = 
9.75/i _1 Gyr, consistent with the lower bound t gc > 10 Gyr for the entire observed 
range of Ho, and yielding a healthy, to > 13 Gyr-old universe for h < 0.75. This 
lower value of Qq is consistent with the dynamical mass estimates from clusters, but 
somewhat below the recent estimates from large-scale flows. Since there is currently 
little theoretical guidance as to why |Ao| is as small as it is, and no firm proof that it 
should vanish, it is probably best to keep an open mind, although the fact that we 
would be living just at the epoch when flo is comparable to Ao might seem to beg for 
explanation. The third possibility is that theoretical prejudice is wrong, and that 
we live in an open, low-density, perhaps purely baryonic universe with negligible 
A, in which case the globular cluster age range is also compatible with somewhat 
larger values of Hq. The challenge in this case is to explain the large-scale flows 
(or attribute them to systematic distance errors) and to form large-scale structure 
without violating CMBR anisotropy constraints. In this case, it also would be an 
odd coincidence that we live just at the epoch when the curvature term in (3.2) is 
becoming appreciable compared to the matter term. 

6. The Early Universe 

We have seen that the FRW models provide a natural framework for under- 
standing observations of the universe on the largest scales, in particular the Hubble 
law. Aside from the Hubble expansion, two other observational pillars support the 
temple of the standard cosmology: the thermal spectrum of the CMBR and the light 
element abundances. Both direct our attention backward to much earlier epochs in 
the history of the universe: the hot big bang. 

6.1. The Cosmic Microwave Background 

The 2.7 K cosmic microwave background radiation, discovered thirty years ago, 
has since been shown to be remarkably isotropic (to roughly one part in 10 5 on angu- 
lar scales larger than a few arc minutes) and to have an exactly thermal blackbody 
spectrum to within the experimental errors. The spectral measurements cover the 
range of photon wavelengths from roughly 0.05 - 50 cm, while the most impressive 
verification of the Planck spectrum comes from the recent COBE FIRAS observa- 
tions in the range 0.05 - 1 cm. The measured FIRAS flux vs. frequency 47 is shown 
in Fig. 6; it is indistinguishable from a Planck spectrum with a temperature of 
T = 2.726 ± 0.010 K. The error bars on the points are too small to see, so they have 
been multiplied by a factor of 100 in the figure. 

The blackbody nature of the CMBR is indicative of a medium that has relaxed 
to thermal equilibrium. However, there are strong grounds for believing that the 
recent universe is transparent to CMBR photons, that is, that the mean free path 
for microwave photon scattering is much longer than the present Hubble length: (i) 
we observe radio galaxies as point sources out to redshifts of order z ~ 4, which re- 
quires microwave transparency out to this scale; (ii) even postulating a fully ionized 
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Figure 6: CMBR brightness vs. frequency (inverse wavelength) measured by COBE 
FIRAS. 

intergalactic medium with closure density (fli gm ~ 1), the mean free time for Comp- 
ton scattering is longer than the Hubble time back to a redshift z <; 10. Therefore, 
the CMBR is not now in thermodynamic equilibrium and could not have been ther- 
malized at recent cosmological epochs. The standard explanation is that the CMBR 
is a remnant of a much earlier hot, dense phase of the universe when the matter 
and radiation were in equilibrium. Indeed, running the Hubble expansion backward 
in time, this is just what one would expect: as the matter is compressed, it heats 
up, becoming a dense, ionized plasma which efficiently scatters and thermalizes the 
radiation on a timescale much shorter than the expansion timescale. 

This argument suggests that the early universe can be accurately described as a 
dilute, adiabatically expanding gas of particles and radiation in local thermal equi- 
librium, uniquely characterized by the instantaneous temperature T(t). While this 
may reflect conditions over much of the early history of the universe, we know it is 
not the full story, because a gas in thermal equilibrium is featureless and structure- 
less, and the observed universe is not so boring. Although thermal equilibrium is a 
good starting point, the interesting epochs in cosmic history-the ones which leave 
potentially observable signatures or relics, such as the light elements, the baryon 
asymmetry, and particle dark matter-are those in which a particle species i goes out 
of thermal equilibrium, because the rate of interactions keeping it in equilibrium, 
Ti, falls below the expansion rate H(t). This happens at a temperature Tp defined 

by 

G)^ 1 ' (6 ' 1) 

Tp is called the freeze-out temperature: at temperatures T < Tp, the particle's 
comoving number density is fixed, n« ~ a~ 3 . Thus, the study of the early universe 
requires both equilibrium and non-equilibrium thermodynamics (the Boltzmann 
equation). For completeness, it is worth noting that there may also be particle relics, 
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such as magnetic monopoles and axions, which were never in thermal equilibrium. 

It is helpful to have some characteristic numbers for the CMBR. For ultra- 
relativistic particles in kinetic equilibrium at temperature T, the energy density 
is 

_^f_ rpi / 1 bosons , 

pi - 3Q ga x | 7/g fermiong t<°.z) 

and the number density 

" = f ^ >< { 3/4 

where ^ is the number of internal degrees of freedom, and the Riemann zeta function 
C(3) ~ 1.202. For blackbody radiation at T = 2.73 K, the present CMBR energy 
density is p 7 = 4.8 x 10~ 34 gm/cm 3 , and the photon number density is n 7 = 420 
cm~ 3 . For comparison, the present baryon density is n# = ps /m p = 1.12 x 10~ 5 f^/i 2 
cm -3 , and the baryon to photon ratio is 

^- = 2.7x ltfn B h 2 . (6.4) 

In the early universe, baryons and photons are kept in thermal contact by Comp- 
ton scattering, e~7 — >• e~7. As the universe expands and cools, eventually it be- 
comes energetically favorable for the ionized plasma to recombine to form neutral 
hydrogen, e~ + p — ► H + 7. When this happens, the density of free electrons drops 
precipitously, and the Compton scattering rate per photon, 7 7 = n e <Jc'omp, falls be- 
low the expansion rate H; as a result, the photons decouple (freeze out) from the 
matter. In the absence of subsequent reionization of the matter, the photons will 
have travelled freely since that epoch, which is therefore called the surface of last 
scattering. 

It is a useful exercise to estimate when hydrogen recombination and photon 
decoupling took place. Since the recombination rate is initially rapid compared to 
the expansion rate, the baryons and photons maintain ionization equilibrium. The 
densities of the non-relativistic baryons in thermal equilibrium are 

O Icy 

(m[F\ 1 (jii — mA 

Ul = 9l v^T ) exp v — t — ) ' ^'^ 

where % = e,p,H, and are chemical potentials. Defining the ionization fraction, 
x = n p /(n p + rifj), and using chemical equilibrium to relate the chemical potentials, 
Pp + Pe = Ph, we arrive at the Saha equation for the equilibrium ionization fraction, 

l-.«,_4^C(3) (rf^f*) , (6.6) 



where B = m e + m p — ran = 13.6 eV is the binding energy of neutral hydrogen. 
Using (6.4) for r\ and taking fish 2 ~ 0.01 from nucleosynthesis, one finds that 
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x eq drops rapidly away from unity-that is, the hydrogen effectively recombines- 
at a temperature T rec ~ 0.3 eV. Since the temperature scales as T = 2.73(1 + z) 
K, this corresponds to the epoch z rec ~ 1300. This result is just what one expects: 
recombination takes place when there are too few photons left in the high-energy tail 
of the thermal distribution to ionize the hydrogen. This happens at a temperature 
somewhat below the binding energy of hydrogen, because the number of photons per 
baryon is so large (77 <C 1). Once x eq drops appreciably below unity, the Compton 
rate per photon, T 7 = xr)n y <Jcomp soon falls below the expansion rate, and the 
photons decouple at a redshift z dec — 1100. From (2.19), this corresponds to a 
decoupling time t dec = 5.6 X lO 12 /^ 1 sec ~ 180, 000/i -1 yr (for J7 = 1). 

If photons decoupled at such an early epoch, why does the CMBR still have such 
a precisely thermal spectrum? The simple reason is that the expansion preserves 
the blackbody thermal spectrum even in the absence of thermalizing interactions. 
At temperatures T > T dec , the photons are in kinetic equilibrium, with a Bose 
distribution function, 

f(p,t<t dec ) = eME 1 /T) _ 1 , (6.7) 

where E = |p|. After decoupling, the photons are free particles, their energies just 
redshifting with the expansion, E(t) = E(t dec )a(t dec )/a(t). However, since they are 
no longer interacting, the occupation number must be conserved for each momen- 
tum state. As a result, the distribution function maintains the form (6.7), with 
effective temperature given by T(t) = T(^ ec )(a(tj ec ) / a(t)); since this holds for each 
momentum state, the spectrum remains thermal. 

In fact, the issue of the preservation of the thermal spectrum is more subtle than 
this argument suggests, and the spectrum observation by FIRAS a more powerful 
tool 48 , because we have not taken into account the possibility that some physical 
process could distort the spectrum away from the form (6.7). For example, annihi- 
lating or decaying elementary particles or exploding black holes could release energy 
into the photon-baryon plasma that is not completely thermalized by the time of 
photon decoupling. Alternatively, after decoupling, an early generation of stars, 
quasars, or active galactic nuclei might release enough energy to heat up and reion- 
ize the baryons, which would then scatter off the CMBR and distort the spectrum. 
The kind of spectral distortion produced depends on the epoch when such a process 
releases energy into the photon-baryon system. Energy released into the baryons at 
z re i ^ 3 x 10 6 is completely thermalized by Compton scattering and bremmstrahlung 
emission: the photons relax to a new equilibrium distribution at a higher temper- 
ature, leaving no distortion. Energy released in the interval 10 5 ^ z re \ ^ 3 X 10 6 
cannot be completely thermalized, because photon production becomes inefficient; 
in this regime, multiple Compton scattering leads to a Bose-Einstein photon distri- 
bution, 
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where the chemical potential ji is related to the energy density injected, Ap, by 
Ap/pf = 0.7///T. The FIRAS results shown in Fig. 6 can be translated into the 
bound |/i/T| < 3.3 x 10~ 4 , severely restricting the energy density injected at such 
times. If energy is released into the electrons later, at z<i ec = 10 3 < z re i < 10 5 , e~j 
scattering produces a Compton y-distortion which can be expressed as a frequency- 
dependent temperature, 

&M ~,J ^ VlT + l ,) , [ -28/ for,/T«l 
T ~ y \Te^l T -l J \yu/T fori//T>l ' K } 

The y-parameter is related to the energy release by y = Ap/4p 7 . Eqn.(6.9) states 
that Compton scattering by hot electrons at T e > T 7 moves photons from low to 
high frequency, depleting the spectrum in the low-frequency Rayleigh- Jeans region 
and enhancing it in the Wien part of the spectrum. Comparison with the data in 
Fig. 6 yields the observational constraint y < 2.5 x 10~ 5 . These bounds together 
severely constrain any model in which energy is released between roughly one and 
300,000 yr after the big bang. 



6.2. The Radiation Era 



Summing eqn.(6.2) over all effectively massless species in equilibrium, we can 
define the total radiation energy density by 

Pr = ^9*T 4 , (6.10) 

where g*(T) counts the number of effective relativistic Bose (B) and Fermi (F) 
degrees of freedom, which may have different effective temperatures, 

9*( T ) = 9i (^) +1^2 9i 

i=B U J 8 i=F 

As the universe expands and cools, the function g*(T) decreases whenever T drops 
below the mass m; of a particle species. For the photon and assuming 3 massless 
neutrinos, it is currently g*(Ta) = 3.36; above the e + e~ mass threshold, g*(l — 
lOOMeV) = 10.75; and at temperatures above the masses of all particles in the 
standard electroweak model, g*(T > 300GeV) <; 100. 

From (6.11), and assuming there are no exotic heretofore unknown relativistic 
particles, the present radiation energy density is p r (t ) = 8.1 x 10~ 34 gm/ cm . 
Comparing with the critical density in eqn.(4.1), the radiation makes a negligible 
contribution to the present density, 

n r = 4.3 x 10~ 5 /T 2 , (6.12) 

confirming our assumption in Sec. 2 that the present universe is matter-dominated. 
On the other hand, since p r .j p m ~ a -1 , the universe was radiation-dominated at 
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early times. From (6.12), the transition from radiation to matter-domination (p r 
p r ) occurred when 



1 + Zeq = ^ = = 2.3 x I0 4 n h 2 , (6.13) 

a eq Tq 



which corresponds to a photon temperature T eq = 5.5f2o/z 2 eV. For Qq = 1, this 
happened when the age of the universe was t eq = 4.4 x 10 10 /i -4 sec. Note that for 
fto = 1 f^s, the universe becomes matter-dominated before (but not long before) 
recombination. 

During the radiation-dominated (RD) epoch (t < t eq ), the spatial curvature term 
in the Friedmann equation can be neglected, since it scales with a more slowly than 
p m or p r and we know that it is at most comparable to the matter density term 
today. The RD era solution is therefore described by (2.20). Combining (2.20) and 
(6.10), we can relate the temperature and expansion rate in the RD epoch: 

H = l.66g 1 J 2 ^- . (6.14) 
Mpi 

Since H = l/2t from (2.20), we can also relate temperature and age: 

0.3M P/ 2.4 ( T \~ 2 

sec . (6.15J 
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Table 2: Thermal history of the universe. 



T(GeV) 


t(sec) 


Event 


10 19 


io- 43 


Planck era: quantum gravity 


10 15 


1Q -35 


Grand unification: inflation, topological defects 


10 3 


io- 12 


Supersymmetry, technicolor 


10 2 


io- 10 


Electroweak transition: SU(2) x U(l) U(l) em 


0.2 


1Q -4.5 


Quark-hadron transition; chiral symmetry breaking 


10~ 3 


1 


Neutrino decoupling 


5 x 10" 4 


4 


e + e~ annihilation 


7 x 10" 5 


200 


Big Bang nucleosynthesis 


10" 9 


5 x 10 11 


Matter-radiation equality 


3 x 10" 10 


10 13 


H recombination; photon decoupling 


6 x 10" 13 


10 17 


Non-linear structures form 


3 x 10" 13 


4 x 10 17 


Chicago Bulls three-peat 






(note: this was given as a prediction at TASI '93) 



A timeline for the early universe with some important epochs delineated is shown 
in Table 2. The CMBR probes conditions back to the time of photon decoupling, 
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while the light element abundances probe conditions at temperatures comparable 
to nuclear binding energies, T ~ MeV. Going backward in time, QCD predicts that 
chiral symmetry should be broken at a temperature of order 100 MeV. At about that 
time, quarks should also become confined inside hadrons, but it is not clear if this is 
a smooth (second-order) transition or a first-order transition involving the release of 
latent heat and the nucleation of hadron bubbles. If the quark-hadron transition is 
first-order, the resulting inhomogeneous baryon distribution might leave observable 
signatures in the light element abundances. 

Going to earlier times, the electroweak symmetry should be restored at a tem- 
perature of order 100 GeV. The dynamics of this transition is currently an area of 
active investigation: if it is a first-order transition, it is possible that the resulting 
non-equilibrium conditions were ripe for the generation of the baryon asymmetry. 
If baryons and anti-baryons had been present in exactly equal numbers in the early 
universe, their annihilation would eventually have driven the baryon to photon ra- 
tio down to a value many orders of magnitude smaller than (6.4). Thus, when the 
baryons were relativistic, there must have been a small asymmetry between the 
density of baryons and anti-baryons, nfi/nj = (n& — ng)/n& ~ r\. As first pointed out 
by Sakharov, the generation of such a baryon asymmetry requires baryon-number 
and CP-violating interactions as well as a departure from thermal equilibrium (since 
rib = ni in equilibrium). It is currently thought that baryogenesis takes place at the 
GUT or electroweak eras. 

Going back earlier than 10~ 10 sec, we must invoke physics beyond the stan- 
dard electroweak model, and the events become increasingly speculative. Particle 
physics models suggest there may be new physics lurking at the TeV scale-perhaps 
supersymmetry, technicolor, or various extensions of the standard model. In the 
simplest grand unified theories, the strong and electroweak interactions are unified 
(the symmetry between them restored) at an energy scale of order 10 15 GeV. A cos- 
mological phase transition at that epoch might lead to inflation or to the generation 
of topological defects such as monopoles, cosmic strings, or textures; the resulting 
density and gravitational wave perturbations produced could provide the seeds for 
large-scale structure and leave a signature in the CMBR anisotropy. Classical cos- 
mology runs into a wall at the Planck era, t ~ 10~ 43 sec: at that epoch, quantum 
fluctuations in the spacetime metric are expected to be large, and a quantum theory 
of gravity is required. If superstrings provide the fundamental description of nature, 
inherently stringy effects would become important around that scale. 

6.3. Relic Neutrinos and Hot Dark Matter 

As an interesting example of interactions freezing out and leaving a relic species, 
consider light neutrinos. The cross-section for v e e~ — >• v e e~ scattering through W 
and Z exchange is of order a v ~ a^p 2 / (p 2 — M 2 ) 2 , where a w is the weak fine structure 
constant, M ~ Mz, Mw is the gauge boson mass, and p are the fermion momenta in 
the center of mass. First consider the very early universe, at temperatures T 3> M. 
In this case, the momenta p ~ T 3> M, resulting in a cross-section (av) ~ o? w jT 2 . 
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Since the particle density n ~ T 3 , the interaction rate is ~ c^T- Comparing 
with the expansion rate, i? ~ T 2 /Mpj, we find that the interaction rate is larger, 
F w /H > 1, at temperatures T ^ o? w Mp\ ~ 10 16 GeV. If we begin our consideration 
of cosmology at the Planck era, T - 10 19 GeV, this tells us that gauge interactions 
are initially out of equilibrium, but local thermal equilibrium is a good working 
assumption below about 10 16 GeV. As a result, inferences about physics near the 
Planck scale which rely on equilibrium thermodynamics should be viewed with 
skepticism. 

The situation changes when the temperature falls below the gauge boson masses, 
T <C M. In this regime, we have M ^> T ~ p 3> m e , and the ve~ cross-section 
becomes (av) ~ a 2 w T 2 /M 4 ~ G 2 F T 2 , where G F is the Fermi constant, G F ~ 1CT 5 
GeV -2 . Comparison with the expansion rate now gives T w /H ~ GpMpjT 3 , which 
falls below unity at the neutrino freeze-out temperature, Tp ~ 1 MeV. For light 
neutrinos, we have Tp 3> m v , so the neutrinos are ultrarelativistic, with a density 
comparable to photons, when they freeze out; such particles are called hot relics, as 
opposed to cold relics, which are non-relativistic (Tp ^ m) when they freeze out. 

Consider the implications of neutrino decoupling for the present neutrino density. 
We will make use of the fact that, in thermal equilibrium, the total entropy in a 
comoving volume is conserved by the first law of thermodynamics, S = so? = 
constant, where the entropy density is 

(6.16) 

Soon after the neutrinos decouple, the electrons and positrons annihilate, e + e~ — >• 
77, at T ~ m e ~ 0.5 MeV, converting their entropy into photons. Since this takes 
place after v freeze out, the neutrinos do not partake of the entropy boost. Before 
annihilation, the e + e~ entropy is s e = 4x (7/8) x (27r 2 /45)T 3 , and the photon entropy 
s 7 = 2 X (27r 2 /45)T 3 . After annihilation, the photon entropy has been increased by 
the factor 



_ 2 + 7/2 _ 11 
s 7 ~ ' 2 " ~ 4 



(6.17) 




which means the photons are heated with respect to the neutrinos, 

(6.18) 

Using the present CMBR temperature from COBE, this implies a present neutrino 
temperature of T v = 1.95 K, and a neutrino density (for a single neutrino species) 

n v = (3/4)(T^/T 7 ) 3 n 7 = (3/ll)n 7 = 115 cm" 3 . (6.19) 

We can also use this to calculate the present energy density in radiation: assuming 
all three neutrino species are massless, we have 

7/4 \ 4/3 

<fc(T ) = 2 + 3x2x-^— J = 3.36 , (6.20) 
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which leads directly to (6.12). 

It is possible that one or more of the neutrinos is massive, the current experi- 
mental neutrino mass limits being m Ue < 9 eV, m v < 270 keV, m Vr < 35 MeV. If a 
neutrino is massive and stable (or long-lived compared to the age of the universe), 
it could contribute more substantially to the density of the universe. From (6.19) 
and (4.1), one finds 

tt v h 2 = — V —H- = — eV . (6.21 

p c h~ 2 92 v ' 

From sec. 5, neutrinos could provide the dark matter and close the universe, fi„ = 1, 
if h ~ 0.5, which would require a mass m v ~ 20 eV. While larger than the v e mass 
limit, either the muon and tau neutrinos could have masses in this range; moreover, 
even an electron neutrino at the upper end of the experimental mass limit would 
contribute substantially to the energy density. 

If neutrinos are the dark matter, an interesting argument due to Tremaine and 
Gunn 43 shows that they cannot cluster in galaxies, so they would not provide the 
dark halos of the kind shown in Fig. 5. To see the point, suppose neutrinos were 
clustered in the halo of a galaxy like the Milky Way. Their number density would 
be n v = J d 3 pf u (p) ~ fu(p 3 ) ~ fu( m u v h) 3 , where the halo velocity dispersion is 
Vh ~ 200 km/sec. By the Pauli exclusion principle, the neutrino occupation number 
must satisfy f u (p) < 1, which implies the constraint n v j{m v v^f' = p v /( m t v h) ^ 1- 
The local halo density in the model of Fig. 5 is of order p/> ~ 0.4 GeV cm -3 ; 
for this to be composed of neutrinos would require m v > (p h /v 3 h ) 1/4 ~ 50 eV. 
Extending this argument to dwarf galaxies leads to the tighter constraint m„ ^ 
(100 - several hundred) eV. These numbers are significantly higher than the closure 
density neutrino mass; turning the argument around, we conclude that in a neutrino- 
dominated universe, galaxy halos must be baryonic. 



6.4. Relic WIMPS and Cold Dark Matter 



In the discussion above, we assumed that the neutrinos are light compared to 
their freeze-out temperature, m v < Tj? ~ 1 MeV. If neutrinos, or some other 
weakly interacting particle, were much heavier, the situation changes: the particle 
becomes non-relativistic before freeze-out, so its abundance relative to radiation 
is depleted by annihilations-such particles are collectively known as cold relics or 
WIMPs. Examples of cold relics are supersymmetric neutralinos, fermionic partners 
of standard model bosons, with masses in the range m ~ 10 GeV —1 TeV. If PL- 
parity is conserved, the lightest supersymmetric particle (LSP) is stable and is a 
prime candidate for (cold) dark matter with flisp ~ 1- 

Accurate calculations of relic WIMP abundances generally require numerical 
solution of the Boltzmann equation, which describes how the WIMP abundance 
relative to its equilibrium value evolves over time in the expanding universe. How- 
ever, one can get an order of magnitude estimate from eqn.(6.1), which states that 
the relic particle abundance (relative to the entropy density) is approximately the 
equilibrium abundance at freeze-out. Using the annihilation rate at equilibrium, 
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r = n eq {<jv), with n eq given by the Boltzmann form (6.5) (here assuming /i = 0) 
and the expansion rate in (6.14), the freeze-out condition becomes (dropping factors 
of 2tt) 



m 

(6.23) 



Thus, the relic WIMP to entropy ratio is 

n eq \ {mT F fl 2 e- m l TF 



t f T F Mpim(av) T F 

For weakly interacting particles with masses in the GeV - TeV range, the ratio 
m/Tp ~ 10 — 20 at freeze-out. For concreteness, consider a massive neutrino 
(other WIMPs are qualitatively similar but differ in quantitative detail): for non- 
relativistic neutrinos of mass m, the four-momentum p ~ m, and the annihilation 
cross-section estimated at the beginning of the last section becomes o v ~ G F m . 
Substitution into (6.23) gives n/s ~ 10~ 8 (m/GeV)~ 3 . Using eqns.(6.10), (6.12), and 
(6.16), this corresponds to a neutrino mass density of order fl„h 2 ~ (m/GeV)~ 2 . 
This tells us that a particle with mass in the GeV range and annihilation cross- 
section typical of the weak interactions will have a relic density comparable to the 
closure density of the universe. Note that such cold particles have no trouble fitting 
into galaxy halos, unlike light neutrinos. 

7. The Inflationary Universe 



7.1. Motivation: the Horizon and Flatness problems 

The inflationary scenario originally arose out of the attempt to solve several 
puzzles which arise when one extrapolates the standard cosmology back to the very 
early universe. The most important of these puzzles are the horizon and flatness 
problems, which we consider in turn. 

Consider a light signal emitted from the origin at time t e ; from eqn.(2.10), it 
will reach a coordinate distance at a time t given by 

r h dr f f dt . 

(7.1) 



Vl - kr 2 Jt e a(t) 

If the second integral converges as t e — ► 0, then there exists a particle horizon: two 
fundamental observers separated by a coordinate distance larger than r/,(t e = 0) 
will not yet be in causal contact at time t. The proper distance from the origin to 
the coordinate r^(t e = 0) at time t is called the particle horizon radius d^, 

d h (t) = a(t) fi . (7.2) 
Jo a(t) 
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In the standard cosmology, in the early radiation- dominated phase, a(t) oc i 1 / 2 , and 
the particle horizon dh ~ ct. Note that this is also comparable to the Hubble length, 
~ t. More generally, dh is finite if the scale factor a(t) grows more slowly 
than t at early times; from (2.21), the criterion for this is p > — p/3, a reasonable 
condition for most fluids. Now recall that, in the absence of early reionization, the 
CMBR photons last scattered at a time t^ ec ~ 5.6 X lO 12 /^ 1 sec. If Qq ~ 1, photons 
emitted from two points separated by the particle horizon distance at the time of 
last scattering would arrive with an angular separation of order 1 degree at the 
observer today. Therefore, in the standard cosmology, when we look at a map of 
the microwave sky, we would expect it to look very anisotropic on angular scales 
larger than a degree, because we are seeing photons emitted from regions which were 
not yet in causal contact with each other at the time of last scattering. Yet the 
microwave sky is remarkably isotropic over all scales larger than a few arcminutes, 
to of order one part in 10 5 . This is the horizon problem: the CMBR is isotropic 
over scales which were not yet in causal contact when the photons last scattered. 
This argument might appear circular, for we have done this calculation using the 
homogeneous and isotropic FRW model, so by assumption the CMBR must appear 
isotropic over all scales, and one might conclude that there is no problem. This 
objection is spurious: we could easily repeat the derivation of the particle horizon 
in a cosmological model which is nearly but not exactly homogeneous and isotropic 
(say, with temperature fluctuations which are of order 10~ 2 instead of 10~ 5 ) and 
end up with essentially the same answer. That is, we do not need to assume exact 
homogeneity and isotropy for there to be a horizon problem. 

The second puzzle of the standard cosmology is the flatness problem. In essence, 
this is the puzzle of why the spatial curvature term in (3.2) does not presently 
dominate over the matter-density term by a large margin: if k ^ and Ao = 0, we 
can rewrite (3.2) as 

WTT\ = almHt) ' (7 ' 3) 

From (2.21), a 2 H 2 = a 2 is a decreasing function of time if w > —1/3. Consequently, 
in the standard cosmology (with non-negative fluid pressure), the observational fact 
that f^o is within an order of magnitude of unity means that Cl(t) must have been 
extremely close to unity at earlier times: at the Planck time, for example, it requires 
|Sl(t_P/) — 1| ^ 10~ 60 . In the absence of other prejudices, one might have guessed that 
the universe would emerge from the quantum gravity Planck era with the matter- 
density and curvature terms having comparable amplitudes, but the result above 
implies that they must have been identical at that time to one part in 10 60 . Another 
way to state this is: if conditions had been just slightly different at the Planck time, 
the universe would have gone into free expansion (a ~ t) or have recollapsed long 
before the present epoch. In other words, the natural timescale for the universe to 
become curvature-dominated is the Planck time, tpi ~ 10~ 43 sec, yet our universe 
is still not strongly curvature-dominated at to ~ 10 17 sec ~ 10 60 tp/. 

It is important to emphasize that these two puzzles are not inconsistencies of 



32 



the standard cosmological model. Rather they point to features of the observed 
universe which the standard model does not explain and which moreover appear 
highly unlikely when the standard cosmology is embedded in a somewhat larger 
class of cosmological models. That is, they are both problems of initial conditions: 
in the "phase space" of initial conditions, the set of initial data that evolve to a 
Universe like ours is very small. A more physical way to say this is that homogeneity 
and spatial flatness are unstable properties of the model, so the present nearly 
homogeneous and spatially flat state of the Universe appears to be very sensitive to 
the initial conditions. Inflation was designed to reduce this sensitivity by widening 
the class of initial conditions which evolve to a nearly homogeneous and spatially 
flat state within the observable universe. This is the sense in which inflation is said 
to 'solve' the horizon and flatness problems. 

7.2. Inflation: necessary ingredients 

The discussion of the horizon and flatness problems implicitly pointed to their 
solution: what is wanted is an early era when the scale factor accelerates, a > 0, for 
this implies that a is increasing with time, i.e., that a grows faster than t. Because 
of the rapid growth of the scale factor when it is accelerating, this was dubbed 
'inflation' by Guth 49 . From (2.13), acceleration requires that the energy density be 
dominated by some fluid with equation of state p < — p/3. How long must such an 
inflationary epoch last? Suppose the accelerated expansion begins at an epoch t; 
and ends at t e . Then, from (7.3), we can compare the present and initial curvature 
terms, 
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In the last expression, I have assumed for the sake of argument that the universe is 
radiation-dominated for t > t e , that the inflation epoch ends around the GUT scale, 
T e ~ 10 15 GeV, and that H e ~ Hi. Here, JV e is the number of e-folds of growth of 
the scale factor during inflation, 



We would like the ratio on the left hand side of (7.4) to be of order unity, so that fl{ 
does not have to be extremely close to one. The right hand side of (7.4) then implies 
that the scale factor must grow by at least N e ~ 60 e-folds during the inflationary 
epoch. If fli ~ 1, this is identical to the condition that our observable universe have 
emerged from a region that was causally connected at the onset of inflation-thus 
solving the horizon problem as well. 

The pleasing feature of inflation is that it drives the spatial curvature to zero 
via rapid expansion of the spatial hypersurfaces. After inflation, the curvature term 
again grows in importance, but it has been reset by inflation to such a tiny value 
that it is still catching up to the matter term today. In fact, if the universe enters an 




(7.5) 
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extended inflationary period, it is quite plausible that it lasts considerably longer 
than the required 60 e-folds of expansion, in which case the spatial curvature term 
would still be exponentially small today In other words, unless JV e is almost exactly 
60, a sufficiently long period of inflation implies that the present universe should be 
observationally indistinguishable from being spatially flat, i.e., J7q + Aq = 1. 

There is a second necessary ingredient for inflation to be a viable scenario for 
the very early universe. The accelerated expansion not only drives the spatial 
curvature to zero, but it also drives the density of matter and radiation to zero. 
At the end of N e > 60 e-folds of expansion, the radiation temperature would be at 
most T e = Tie~ 60 which is less than 100 eV if T- < Mpi. This would invalidate the 
argument in eqn.(7.4), where we assumed T e ~ Mqut, an d, among other things, 
make it impossible to account for the light element abundances through primordial 
nucleosynthesis. Thus, a period of accelerated expansion by itself is not enough to 
solve the horizon and flatness problems: at the end of inflation, the energy density in 
the fluid that drives the inflationary expansion must be converted into radiation, to 
repopulate and reheat the universe. A minimal condition for successful inflation is 
that the reheat temperature after inflation be high enough for the baryon asymmetry 
to be generated; if baryogenesis takes place at the electroweak scale, this implies 
T e > 100 GeV. 



7.3. Inflation: Scalar Field Dynamics 

In 1980, Guth 49 proposed that a scalar field trapped in a 'false vacuum' state 
with non-zero potential energy could act as the driving force of inflation and thereby 
solve the horizon and flatness problems. It is easy to see that this fits the bill for 
accelerated expansion: the energy-momentum components for a classical scalar field 
0(x, t) with potential V(4>) are 

P4, = T 6 ^ ( ) 

and 

If the field is trapped in a local minimum 4> c of the potential, it will relax to a static 
configuration with p$ = —p^ = — V(<f> c ) = constant, which satisfies the criterion 
for accelerated expansion. In this case, eqn.(2.12) becomes H 2 = 87rGV((j) c )/3 
= constant, and the solution is the exponential de Sitter expansion of (3.7). In 
other words, a field dominated by constant potential energy acts as an effective 
cosmological constant. 

While a promising idea, 'false vacuum' inflation did not satisfactorily incorporate 
the second necessary ingredient of reheating: once the field is trapped for sufficiently 
long in the false vacuum state, it must tunnel through the potential barrier which 
held it there, in order to reach the true vacuum where V(<f>t) = 0. In this case, 
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all the potential energy of the false vacuum ends up in the walls of the nucleated 
bubbles of true vacuum, and is not converted into radiation by bubble collisions 
efficiently enough. The remedy was soon provided by Linde 50 and by Albrecht 
and Steinhardt 51 : instead of the field being trapped in a metastable minimum of 
the potential, it can be classically evolving on a very flat potential if it is initially 
displaced from the potential minimum. If the potential is sufficiently flat, the kinetic 
and potential energy terms in eqns. (7.6-7) will redshift away with the expansion, 
leaving the potential term to dominate as before. Inflation ends when the field 
reaches a steeper part of the potential, the field speeds up and eventually oscillates 
about its potential minimum. If <j) is coupled to lighter particles, these coherent field 
oscillations lead to particle creation and thereby reheat the universe: the scalar field 
energy is converted to radiation. This 'new inflation' model thus solves the 'graceful 
exit' problem of the original false vacuum scenario. 

In these early inflation models, it was assumed that the scalar field driving 
inflation, the inflaton, was a Higgs field associated with the spontaneous breakdown 
of the grand unified gauge symmetry. The field would be initially displaced from 
its global minimum by finite temperature effects. It was soon realized, however, 
that the inflaton must be extremely weakly self-coupled in order for its quantum 
fluctuations to generate an acceptable amplitude of density perturbations. Since 
Higgs fields are coupled to gauge fields, however, radiative corrections typically 
generate large self-couplings. As a result, the concept of inflation was divorced 
from Higgs fields and gauge symmetry breaking, and a pandora's box of particle 
physics models was opened up. Nevertheless, the different models of inflation share 
many common dynamical features, and it is those which I will focus on. 

To see how inflation works in a little more detail, consider the dynamics of the 
evolving scalar field. Within each Hubble volume, (i.e., ignoring spatial gradients) 
the evolution of the field is described by the classical equation of motion for a 
homogeneous field <f>(t), 

4> + 3H<f> + T<f> + V'(<f>) = , (7.8) 

where T is the decay width of the inflaton field, a phenomenological term which has 
been introduced to roughly model the back-reaction effects of particle creation and 
reheating on the scalar field. The expansion rate H = a/a is determined by the 
Friedmann equation, 

8tt 
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v(4>) + 



(7.9) 
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and it is also useful to have the second order Friedmann equation, 



From (7.10), the condition that the universe be inflating, a > 0, requires that the 
kinetic energy by sub-dominant, <f> 2 < V. A sufficient, although not necessary, 
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condition for this is that the field be slowly rolling (SR) in its potential. In this 
context, 'slowly rolling' is a technical term: the field is said to be slowly rolling 
when its motion is overdamped, i.e., 

4><.$H4> , (7.11) 

so that the <f> term can be dropped in the equation of motion (7.8) (we assume 
r <C H during this phase). From (7.8), the defining SR condition implies that 
<fi 2 <C 2V(<f>); thus, if the SR condition is well satisfied, the universe is inflating. 

Slow- rollover is analytically convenient, for it implies two consistency conditions 
on the slope and curvature of the potential: 

£748^ , (7.12) 

where the slow- roll equation of motion has been used in the last equality (eqn.(7.8) 
without the second derivativer term). Suppose inflation begins when the field value 
is <p(ti) = <f>i, and the SR epoch ends when 4> reaches a value <f> e , at which one of 
the inequalities (7.12) is violated. To solve the cosmological puzzles, we demand 
that the scale factor of the universe inflates by at least 60 e-foldings during the SR 
regime, 

t u -Sty r<t>e V(d>) 

N(<f>i,(f> e ) = \n(a e /oi) = Hdt = — ^ / j^±d<f> > 60 . (7.13) 

Jt, Mpi J<p t V\<p) 

Once the potential is chosen, (7.12) determines 4> e in terms of the parameters in V. 
The condition (7.13) then constrains the initial value ^ of the field. 

Once (j) grows beyond <f> e , the SR condition breaks down, and the field evolution is 
more appropriately described in terms of oscillations about the potential minimum. 
These coherent oscillations excite the creation of particles to which <f> is coupled, 
which in turn damps the oscillations and reheats the universe once the created 
particles thermalize by scattering. By energy conservation, the phenomenological 
damping term T in the scalar equation of motion-the scalar decay rate-acts as a 
source for radiation, 

p r . + mp r = T^ 2 . (7.14) 

We can identify two qualitative regimes for reheating: (i) if T > H, reheating is 
very efficient: nearly all the scalar potential energy is converted to particles, leading 
to a reheat temperature given by T^ H ~ V((f) e ); (ii) on the other hand, if T < H, 
reheating is inefficient, and the scalar field oscillations partially redshift away before 
converting to particles. In this case, the reheating temperature is lower, and is given 
by 

Trh = (45/47r 3 ^) 1 ^ 4 yT^Mp/ , (7.15) 

where is the number of relativistic degrees of freedom. Clearly the reheating 
efficiency is determined by the couplings of the inflaton to other fields. Efficient 
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reheating, say, with a reheat temperature high enough to allow baryogenesis at 
the GUT scale, Tjhj ~ 10 14 GeV, requires relatively strong couplings. In many 
models, however, such couplings induce loop corrections to the scalar self-interaction 
which would upset the requisite flatness of its potential. As a result, a high reheat 
temperature is not always possible. 

7.4. Quantum Fluctuations and Density Perturbations 

Soon after the arrival of false vacuum inflation, it was appreciated that inflation 
could in principle 52 ' 53 provide one of the holy grails of cosmology: a causal mecha- 
nism for the origin of density fluctuations that later grow to form large-scale struc- 
ture. The problem can be seen by first considering the standard cosmology without 
inflation: in this case, the physical wavelength associated with a density perturba- 
tion mode grows with the scale factor, X p h ys ~ a (t), while the Hubble length scales 
up more rapidly, H _1 ~ t (recall the current Hubble distance is ~ Hq 1 = 3000/i _1 
Mpc). Consequently, at a redshift larger than z\, where 1 + z\ ~ (Hq 1 /Xq) 2 , per- 
turbations with present scale Aq were larger than the instantaneous Hubble radius, 
which sets the lengthscale for causal processes. For perturbations on cosmologi- 
cally interesting scales today, say, Ao ~ 1 — 3000/i _1 Mpc, this implies z\ ^ 10 6 , 
or T\ ^ 100 eV (here, for the sake of simplicity, I have assumed matter-dominated 
expansion throughout; the reader can insert the appropriate radiation-dominated 
corrections to make the numbers more rigorous). In the standard cosmology, how- 
ever, curvature perturbations cannot be causally generated on scales larger than the 
instantaneous Hubble radius, because of local energy conservation. Thus, pertur- 
bations on scales of galaxies and larger must have been generated by some physical 
process acting at T ^ 100 eV, after these scales crossed inside the Hubble ra- 
dius. While some work has been done on a 'late time phase transition' at such a 
scale, 54 such models do require new physics at surprisingly low energy scales, and 
have not been met with overwhelming enthusiasm. Moreover, the resulting CMBR 
anisotropics in such scenarios may be intolerably large . If one does not postulate a 
late-time origin for large-scale perturbations, then they must simply be postulated 
as initial conditions, because they could not have had a causal origin at epochs 
z ^ z\. (Note that this argument applies to curvature perturbations; an important 
loophole is the early generation of large-scale isocurvature perturbations.) 

An early epoch of inflation alters this situation: during accelerated expansion, 
the Hubble radius H~ x grows more slowly than the scale factor (in a de Sitter epoch, 
H^ 1 is constant). Thus, a perturbation of comoving wavelength A could be created 
causally on a physical scale less than during inflation, expand outside the 

Hubble radius at some time ^(A) during inflation, and then eventually 're-enter' 
the Hubble radius in the recent radiation- or matter-dominated era at a time ig(A). 

This hypothetical possibility was turned into a physical plausibility when it was 
realized that quantum fluctuations of the slowly rolling field in new inflation could 
provide a mechanism for generating density perturbations 56 . In a nutshell, one 
treats the spatial fluctuations of the field about its homogeneous classical expec- 
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tation value as an approximately massless quantum field in de Sitter space; in its 
vacuum state, such a field has a spectrum of zero-point fluctuations, with an rms 
amplitude A(f> = H j2n. These spatial variations in <f> correspond to fluctuations in 
the inflaton energy density, which are converted into adiabatic density fluctuations 
in all species during reheating. The resulting density perturbation amplitude when 
a comoving wavelength A re-enters the Hubble radius is given approximately by 



It is straightforward to show that scales corresponding to the range from the galaxy 
scale to the present Hubble radius, Ao = 1 — 3000/i _1 Mpc, crossed outside the Hubble 
radius in the interval 50 - 60 e-folds before the end of inflation. During slow- rollover, 
the inflaton is effectively at its terminal velocity, which changes slowly over time; 
since the period of ten e-folds of the scale factor is brief, the right hand side of 
(7.16) changes very little over this interval. As a result, the density perturbation 
amplitude at Hubble radius crossing (t#(A)) is nearly constant. 

Now, the density perturbation amplitude at Hubble radius crossing is essentially 
the gravitational potential perturbation on that scale, S&\ ~ GSM\/\ ~ G5p\\ 2 ~ 
(Sp/ p)x(Xph ys / H~ v ) 2 , where we have used the Friedmann equation to relate Gp ~ 
H 2 . Thus, inflation predicts that the gravitational potential fluctuations are nearly 
independent of scale, 5<&\ ~ constant. Such a scale-invariant spectrum of potential 
fluctuations was first proposed on different grounds by Harrison, Zel'dovich, and 
Peebles and Yu, ten years before inflation. Note that the gravitational potential 
also sets the scale for the large-angle CMBR anisotropy through the Sachs- Wolfe 
effect: photons climbing out of deeper potential wells at the time of last scattering 
suffer a larger gravitational redshift, resulting in the temperature shift ST/T = 
~ (Sp/ p)\ =H -i (this result holds for the Einstein-de Sitter model, J7q = 1). 
Considerable excitement was generated when analysis of maps of the CMBR sky 
from the first year of COBE DMR data yielded a large-angle anisotropy consistent 
with the scale-invariant spectrum 3 . 

To get a feel for the result (7.16), I sketch out the steps of the calculation here 
(for a particularly complete discussion, see 57 ). The idea is to treat small fluctuations 
semi- classically, by expanding the scalar field and metric perturbatively around their 
homogeneous solutions, e.g., 



where (f> (t) is the classical homogeneous solution to the scalar equation of motion in 
the FRW background, eqn.(7.8), and 8<f> is a quantum field operator with equal-time 
commutation relation 



A similar expansion is performed for the metric. It is convenient to consider the 






(7.17) 



[6<t>(x,t),d6<t>(x!,t)/dt] = ^£ 3 (x-x')/a 3 (^) 



(7.18) 
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Fourier transform 



X ' t} = / jBf K^(t)e 4k ' x + h.c] , (7.19) 

where a^ is the annihilation operator for mode k, ipk{t) is the classical mode function, 
and k = ak p h ys is the (fixed) comoving wavenumber of the mode. We will consider 
the power spectrum V<f,(k) of the field, the contribution per logarithmic wavenumber 
interval to the variance of the fluctuation, 

m?) = j T V,(k) . (7.20) 

The power spectrum is proportional to the Fourier transform of the two-point cor- 
relation function of the field in the vacuum, 

V*{k) = d 3 ^ k ' x (0|^(x)^(0)|0) = |^|^| 2 . (7.21) 

In general the mode equation for ip k can be complicated, since it involves coupled 
scalar and gravitational degrees of freedom. However, at least in the synchronous 
gauge (<5goo = 0) ; the term involving metric perturbations can be neglected, yielding 

k 2 

iP k + 3H ip k + —ip k + V"(<P o )iP k = . (7.22) 
er 

The slow rollover conditions imply that the potential term in (7.22) is also negligible, 
reducing the equation of motion to that of a free massless field in de Sitter space. 
The normalized solution is 

Jka for . (7.23) 

constant for k/aH <C 1 

The mode crosses outside the Hubble radius during inflation when X p h ys = H^ 1 , 
i.e., when k = aH. The solution (7.23) states that the mode oscillates inside the 
Hubble radius and then freezes out when it crosses outside. This is the mathematical 
statement of the physical intuition that microphysics only operates coherently on 
scales less than Hubble radius. For modes inside the Hubble radius, from (7.21) 
and (7.23) the power spectrum can be written V^{k) = (H /2Ti) 2 (k / Ha) 2 . When a 
mode crosses outside the Hubble radius, the rms fluctuation on that scale is given 

by 

(A0) M = Vl /2 (k = aH) = ^ . (7.24) 

To find the resulting density perturbation amplitude, one must solve the per- 
turbed Einstein equations; due to the gauge freedom in the choice of time coordi- 
nate, that is, in the way spacetime can be sliced up into spacelike hyp ersurf aces, one 
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chooses a gauge to simplify the problem at hand. It is particularly convenient to 
make use of a gauge-invariant measure of the perturbation amplitude for adiabatic 
(curvature) perturbations, 

C = 3$ + ^- , (7.25) 

Po + Po 

where $ is an analog of the Newtonian potential that completely specifies the in- 
trinsic Ricci curvature of the perturbed 3-surfaces. The variable £ is useful because 
it generally satisfies ( = constant for perturbations far outside the Hubble radius, 
k <C aH, both during and after inflation. Moreover, the potential $ is negligible 
at Hubble radius crossing, so, defining the power spectrum of £ by analogy with 
(7.20), we have 

~ Vl'\k) = (Wj-sZiM = , (7 . 26 ) 

t B w ^ \^4>) tA {x) 

where we have used (7.24) and the slow-roll equation of motion, 3H(f> = —V'(<f>). 
This completes the sketch of the derivation of (7.16). 

Using the slow-roll equation of motion, the perturbation amplitude can be ex- 
pressed in terms of the height and slope of the potential, 

T/3/2 

^ = ^vW, ■ <7 - 27) 

A useful way to characterize the scale-dependence of the perturbations is to consider 
the density perturbation power spectrum, P p (k) = \5k\ 2 , where 8j. is the Fourier 
transform of <5p(x)/p, and we relate different conventions for what is called the 
'power spectrum' by using the notation V = k 3 P. The density spectrum can be 
written as a power law in wavenumber, P p {k) ~ k n ", where n s = 1 corresponds 
to the scale-invariant spectrum, i.e., to Vq = constant. (To see this, recall that 
S&\ ~ (Sp/ p)\(\/H~ 1 ) 2 , and that the contribution per logarithmic wavenumber 
interval to the rms fluctuation on scale A is (Sp/p)\ ~ k 3 / 2 \ 5k \ . ) Small deviations 
from scale-invariance can be expressed in terms of the index n s — 1. 

One can carry through an argument similar to that above for gravitational wave 
perturbations in de Sitter space. Expanding the metric about the de Sitter FRW 
solution, the tensor modes /i + x satisfy the massless scalar field equation (7.22) 
(with no potential term). Defining a canonically normalized scalar field for each 
polarization state, h +jX = ^/\Q^lG(f) + ^, the tensor amplitude at Hubble crossing is 

2 H V 1 ! 2 



l k=Ha 
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The induced Sachs- Wolfe anisotropy on large scales is 5T/T ~ hk-jj a . 



(7.28) 
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7.5. Chaotic Inflation: a worked example 



The result (7.26) for the perturbation amplitude provides a severe constraint on 
the form of the inflaton potential. To get a feeling for the numbers, consider one 
of the simplest forms for the potential, a pure quartic V(<f>) = A(/> 4 /4!, an example 
of Linde's chaotic inflation. In this case, the field rolls down to the origin from its 
initial value; the slow roll conditions (7.12) are satisfied for <p > <f) e = Mp//v / 2~7r, 
and the number of inflation e-folds is given by (7.13), 

^•« = ]g;(^) = Sf-i ■ ("9) 

Thus, the condition of sufficient inflation, N e > 60, requires <f>i > 3.1Mp/, and the 
density perturbation amplitude on the scale of the present Hubble radius is 

Pi' 2 = (vrA) 1 / 2 (^-) 3 - 50A 1 / 2 . (7.30) 

The perturbation spectrum is almost exactly scale-invariant in this model, n s ~ 
0.95, because (f> changes very little in the interval between 60 and 50 e-folds before the 

1 /2 A 

end of inflation. The COBE observation of CMBR anisotropy indicates P^ ~ 10 , 
which implies the inflaton self-coupling must satisfy 

A < 4 x 10~ 12 . (7.31) 



7.6. Naturalness and the Small Coupling Problem 

The preceding example illustrates a general conclusion that applies to all in- 
flation models: the inflaton field must be very weakly self-coupled. In different 
models, this constraint appears in different guises, but it is alway present in one 
form or another: the potential must contain a very small dimensionless number of 
order 10~ 12 . 

Attitudes concerning this problem vary widely among inflation theorists: it is 
often said that A^ ~ 10~ 12 constitutes unacceptable 'fine tuning', but that is a 
misapplication of the fine tuning concept in this context. To others, it is not an issue 
of great concern, because we know there exist other small numbers in physics, such 
as lepton and quark Yukawa couplings gy ~ 10~ 5 and the ratio M wea f./Mpi ~ 10~ 17 . 
Partly as a consequence of the latter view, in recent years, it has become customary 
to decouple the inflaton completely from particle physics models, to specify an 
'inflaton sector' with the requisite properties, with little regard for its physical 
origin. 

Nevertheless, it is meaningful and important to ask whether such a small value 
for A^ is in principle unnatural. Clearly, the answer depends on the particle physics 
model within which <fi is embedded and on one's interpretation of naturalness. A 
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small parameter A is said to be "technically natural" if it is protected against large 
radiative corrections by a symmetry, i.e., if setting A — ► increases the symmetry 
of the system. For example, in this way, low energy supersymmetry might protect 
the small ratio M wea k/Mpi by cancelling boson and fermion loops. However, in 
technically natural models, the small coupling A^, while stable against radiative 
corrections, is itself unexplained, and is generally postulated (i.e., put in by hand) 
solely in order to generate successful inflation. Technical naturalness is a useful 
concept for low energy effective Lagrangians, like the electroweak theory and its 
supersymmetric extensions, but it points to a more fundamental level of theory for 
its origin. For example, the underlying origin of the mass hierarchy M wea k/Mpi in 
supersymmetric theories is thought to be associated with hidden sector physics in 
supergravity theories, which start out with no small dimensionless couplings at the 
Planck scale. Since inflation takes place relatively close to the Planck scale, it would 
be preferable to find the inflaton in such particle physics models which are "strongly 
natural", that is, which have no small numbers in the fundamental Lagrangian. 

In a strongly natural gauge theory, all small dimensionless parameters ulti- 
mately arise dynamically, e.g., from renormalization group (or instanton) factors 
like exp( — 1/a), where a is a gauge coupling. In particular, in an asymptotically 
free theory, the scale M\, at which a logarithmically running coupling constant 
becomes unity, is M\ ~ il^e -1 /", where Mi is the fundamental mass scale in the 
theory. In some models, the inflaton coupling A^ arises from a ratio of mass scales, 
A^ ~ ( .\/| As a result, in such models, A^, is naturally exponentially sup- 
pressed, A^ ~ e~ q l a . 

One example of a class of such models is 'natural inflation', 58 ' 59 in which the role 
of the inflaton is played by a pseudo-Nambu-Goldstone-boson (in particle physics 
models, an example of this is the axion). In its simplest incarnation, consider a 
global £7(1) symmetry spontaneously broken at a high-energy scale / by a non- 
zero expectation value for a complex scalar field, the potential for which takes the 
familiar Mexican-hat (or wine bottle) form. At scales much below /, the only 
relevant degree of freedom is the massless angular variable <f> around the bottom of 
the potential. At an energy scale A <C /, the symmetry is explicitly broken-the 
Mexican-hat is tilted, generating a potential for <j) which is generally of the form 

V(<f>) = A 4 [l±cos(iV<£//)] . (7.32) 

In axion models, A characterizes the scale at which a running gauge coupling 
becomes strong, and is related to the fundamental scale / by A ~ fe~ ga - For 
/ ~ Mpi ~ 10 19 GeV and A ~ M GUT ~ 10 15 GeV, this field can drive inflation; in 
this case, the effective quartic coupling is A^ ~ (A//) 4 ~ 10~ 13 , as required. In the 
standard model and in typical grand unified theories, the only strong gauge coupling 
is that of QCD, with Aqcd — 100 MeV. However, in supergravity and superstring 
models, there is a 'hidden' gauge sector which interacts only gravitationally with 
ordinary quarks and leptons. It has been suggested that the hidden sector gauge 
interactions become strong at a scale comparable to that above, A ~ 10 14 GeV, in 
order to break supersymmetry at the TeV scale in the observable sector. 
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An interesting possibility in this case is that the perturbation spectrum can 
deviate significantly from scale-invariance: the spectral index is given by 



= (7-33) 

The constraint that the universe reheat after inflation to a temperature Trh ^100 
GeV leads to the constraint / jMpi > 0.3, which implies n s <: 0.5, allowing a small 
but significant break from scale-invariance (n s = 1). Contrary to some erroneous 
statements in the literature, such a breaking from scale-invariance does not require 
any fine-tuning of the parameters in this model. Arguments from the formation 
of structure, combined with the COBE observations, at any rate indicate that the 
spectral index n s j£ 0.6 — 0.7. 

8. A Brief Look at Structure Formation 

Let us finish with a cursory look at how primordial perturbations, e.g., from 
inflation, end up looking today. As discussed above, the primordial spectrum from 
inflation is generally a power-law in wavenumber, P p (k) = \8 k (ti)\ 2 = Ak n ", and 
the present spectrum is related to the primordial one through a transfer function, 
\8 k (t )\ 2 = T 2 {k)\8 k {U)\ 2 . The transfer function T(k) encodes the scale-dependence 
of the linear (rf < 1) gravitational evolution of the perturbation modes; it depends 
on the nature of the dark matter (hot or cold), its density (Qbm) as well as Qb and 
Hq. On scales which enter the Hubble radius after the universe becomes matter- 
dominated, k <C k eq ~ 0.2Vtoh 2 Mpc -1 , all perturbations undergo the same growth 
rate, so the transfer function T 2 (k) ~ 1; on smaller scales, k k eq , it bends over 
to T 2 (k) ~ k~ 4 as k — ► oo, reflecting the suppressed growth of fluctuations which 
cross inside the Hubble radius while the universe is still radiation-dominated. For 
standard cold dark matter (CDM), we have h = 0.5 and assume negligible baryon 
density, Qb <C 1, Q co id = 1? leading to the characteristic scale k eq ~ 0.05 Mpc -1 . A 
reasonable analytic fit to the linear transfer function for CDM models is given by 60 



1.7k 9k 3 ^ 2 
Oh + (Qhf/ 2 ' (Qhf 



(in the approximation Qb = 0, = Q co id)- For hot dark matter, on the other hand, 
in the absence of seed perturbations such as cosmic strings, the transfer function is 
also exponentially damped by relativistic neutrino free-streaming for wavenumbers 
larger than k v ~ 0.1(m I/ /20eV) Mpc -1 . Finally, the present density spectrum is 
related to the galaxy power spectrum by a bias prescription; for the simplest linear 
bias model, 

PUk) = h 2 gal P p (k) = b 2 gal T 2 (k)\8 k (t i )\ 2 . (8.2) 

As mentioned in Sec. 2, galaxy catalogs such as the CfA, the IRAS 1.2 Jy and 
QDOT, and the angular APM and EDSGC catalogs, suggest that the spectrum 
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Figure 7: The APM galaxy angular correlation function of Fig. 1 is shown with two 
model predictions: standard (flh = 0.5) CDM and a best-fit low-density CDM model 
(flh = 0.17) (from ref. 62). The upper and lower bracketing curves for each model 
indicate the expected 1 — o spread due to cosmic variance. 

predicted in (8.1) and (8.2) with standard CDM (flh = 0.5), linear bias, and n s = 1 
may not have exactly the right shape to match the observed galaxy spectrum, hav- 
ing relatively too little (too much) power on large (small) scales. Solutions to this 
extra-power problem can be classified according to which element on the right hand 
side of Eqn.(8.2) one tinkers with. In all cases, the aim is to flatten the shape of the 
spectrum at intermediate wavenumbers k ~ 0.05 h Mpc -1 by increasing the relative 
power on these scales compared to smaller wavelengths. For example, as (8.1) indi- 
cates, one can abandon 'standard' CDM and increase the characteristic wavelength 
where the transfer function T(k) bends down, providing more relative large-scale 
power, by reducing Qq from unity; in the context of inflation, this requires either 
a special choice of the number of e-folds N e or the introduction of a cosmological 
constant to ensure Qq + Ao = 1. An example of this is shown in Fig. 7, which again 
shows the APM galaxy angular correlation function plotted against the prediction 
for w(9)/(b ga i(Tg) 2 (see below) in linear theory for standard and low-density CDM, 
both with n s = 1. Recall that lowering Qq also increases the expansion age HqIq. 

Another recently popular alternative is to employ a mixture of hot and cold 
dark matter, with flhot — 1 — fl co id = 0.3; from eqn.(6.21), this can be achieved 
with a light neutrino with mass m v ~ 7 eV. Due to free-streaming, the neutrino 
admixture partially suppresses power on small scales, which is what is wanted. A 
third possibility is to retain standard CDM but consider primordial perturbation 
spectra with n s < 1 from inflation, but it appears that such tilting of the spectrum, 
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when normalized to COBE, does not by itself solve the power problem without 
violating other constraints. Alternatively, one can admit a more complex scheme for 
biasing in which the bias factor in (8.2) is scale-dependent, b g — >• b g (k), and increases 
at large scales 61 ' 62 . Unlike the other fixes, in this case the extra power on large scales 
is an optical illusion, a property of the galaxy field but not the underlying density 
field. This difference leads to qualitatively different behavior of the higher-order 
(e.g., three-point) galaxy correlation functions, an effect which should be testable 
against observations 63 . 

Since all inflation models predict a deviation from scale-invariance at some level, 
it is interesting to see what limits on the spectral index n s follow from COBE in the 
context of standard CDM. For these models, the density power spectra form a two- 
parameter family characterized by the spectral index n s and the normalization A. 
Instead of A, it is common to normalize spectra by the rms linear mass fluctuation 
in spheres of radius 8 h~ x Mpc, erg = ((eiM/M) 2 )^ 8/i _ lM c , where 



filters out the contribution from small scales. Redshift surveys of optically selected 
galaxies (in particular the CfA and APM surveys) indicate that the variance in 
galaxy counts on this scale is of order unity. Thus, in a linear bias model (and 
ignoring non-linear gravitational effects), the bias factor for these galaxies would be 
b op t ~ l/""8j f° r other galaxy populations, b ga i<7% may differ from unity. 

From the Sachs- Wolfe effect, for a given normalization and spectral index, we 
can calculate the expected large-angle CMBR anisotropy. The rms temperature 
fluctuation on the scale of 10° observed by COBE, <r T (10°) = 1.085 x 10" 5 (1 ±0.169), 
then yields a relation between them, erg ~ e - 2 63 ( 1 -™») [i -|- 0.2]. Thus, models which 
obtain more relative large-scale power by tilting the spectrum to n s < 1 reduce 
the power on galaxy scales. If erg 0.5, the power on small scales is reduced to 
the point that galaxy formation would occur too recently; in the context of these 
models, this leads to the constraint n s ^ 0.6 — 0.7. This bound assumes that 
gravitational waves make a negligible contribution to the CMBR signal. While 
this is a good approximation for models such as natural inflation, there are other 
inflation potentials for which the gravity wave contribution to the large- angle CMBR 
anisotropy becomes appreciable when n s < 1. An example is power law inflation, 
produced by an exponential potential for 4>: in this case, the relative contribution 
of the tensor modes compared to the scalar (density perturbation) modes to the 
large-angle anisotropy is R ~ 6(1 — n s ). However, contrary to some statements in 
the literature, this relation is not generic to inflation, but only holds approximately 
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and the window function 





45 



for a restricted class of models. Clearly, for models in which it does apply, the effect 
is to lower the COBE constraint on the perturbation amplitude og for fixed n s . For 
models with n s < 1, this makes the epoch of structure formation more recent, and 
the corresponding lower bound on n s even tighter. As this subject has become a 
booming theoretical industry in the last year, I refer the reader to recent issues of 
Physical Review D for more details. 
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Appendix A: Notation and Units 

It is useful to be able to translate back and forth between astronomical units 
and high energy physics units. The fundamental unit of distance in extragalactic 
astronomy is the parsec (pc): 1 pc = 3.26 light-years = 3.09 x 10 18 cm. The typical 
galactic scale is measured in kiloparsec (kpc): for example, the solar system is 
approximately Rq = 8 — 8.5 kpc from the center of the Milky Way. The nearest large 
galaxy to our own is Andromeda, at a distance of roughly 700 kpc, and the typical 
distance between neighboring large galaxies (outside clusters) is a few Megaparsec 
(Mpc). Absolute extragalactic distances are uncertain by a factor of order 2, and 
are usually given in terms of hr x Mpc, where the Hubble parameter Hq = lOO/i 
km/sec/Mpc, and observations indicate 0.4 < h < 1 (see below). Cosmological 
distances are often estimated using redshifts via the Hubble law, v = cz = H^d, and 
so are sometimes expressed in terms of recession velocity: 100 km/sec = 1 hr x Mpc. 
Astronomical masses and luminosities are measured in solar units, M Q = 1.989 x 10 33 
gm, L Q = 3.86 x 10 33 erg/sec. Large galaxies commonly have a luminosity of order 
10 10 — 10 n L Q and mass of order 10 11 — 10 12 M Q . In particle physics, it is often 
convenient to calculate in natural units, in which h = c = k% = 1. In this case, 
all dimensionful quantities can be expressed in units of energy, say in GeV; for 
example, Length = 1/Energy. At the end of a calculation, to convert an expression 
to physical units, insert appropriate powers of he = 2 x 10~ 14 GeV cm, and convert 
lengths to times by dividing by c = 3 x 10 10 cm/sec. When gravity is involved, use 
the Planck mass (energy) Mpi = G _1//2 = 1.2 x 10 19 GeV. Finally, when converting 
energy to temperature (E = ksT) and vice versa, use 1 eV = 1.16 x 10 4 K. To 
go from particle physics to astronomical units, it is sometimes handy to remember 
that the sun contains about 10 57 protons, i.e., 10 57 GeV ~ 1M Q . 
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